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Abstract: We compute the master integrals for massless two- loop vertex graphs with three off- 
shell legs. These master integrals are relevant for the QCD corrections to H — > V*V* (where 
V = W, Z) and for two-loop studies of the triple gluon (and quark-gluon) vertex. We employ the 
differential equation technique to provide series expansions in e for the various master integrals. 
The results are analytic and contain a new class of two-dimensional harmonic poly logarithms. 
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1. Introduction 

It is well known that many of the loop integrals that appear in Feynman diagram calculations can be 
expressed in terms of hypergeometric functions with parameters that depend on the number of space 
time dimensions d and a number of kinematic scales. However, expressing these hypergeometric 
functions as expansions in e = (4 — d)/2 is rather non-trivial. In general the coefficients involve 
polylogarithms, both of the Nielsen [1] and harmonic [2, 3] varieties, and often new polylogarithms 
need to be introduced. 

Integration by parts [4, 5] (IBP) and Lorentz invariance [6] (LI) identities are crucial in reduc- 
ing [7] the number of master integrals (MI) that actually need to be evaluated and several powerful 
tools have been established to deal with the problem. Often, these methods rely on the link between 
the hypergeometric functions that yield (nested) sums and their integral representations that yield 
polylogarithms. Two of the most powerful analytic methods are the Mellin-Barnes technique [8] 
and the differential equations approach [9] . Both have been used extensively to provide expansions 
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in e for two-loop box graphs with massless internal propagators when all the external legs are on- 
shell [10, 11, 12, 13, 14] and when one of the external legs is off-shell [15, 16]. All of these expansions 
have been checked by Binoth and Heinrich's numerical program [17] for evaluating loop integrals. 
Once analytic expansions for the two-loop master integrals were known, they were rapidly exploited 
in the calculation of the amplitudes for physical scattering processes. 

The two-loop helicity amplitudes have been evaluated for the gluon-gluon [18, 19], quark- 
gluon [20, 21] and quark-quark [22] processes and have confirmed the earlier "squared" matrix 
elements [23, 24, 25, 26]. Similarly, amplitudes for the phenomenologically important gg — > 77 [27] 
and qq — ► 77 [28] processes as well as 77 — > 77 [29, 30] and (massless) Bhabha scattering [31] 
have also been computed. The processes with one off-shell leg include e + e~ — > 3 jets [32, 33, 34] 
which is crucial in making a precise determination of the strong coupling constant at the NLC. 
Progress is also being made in calculating the two-loop QED corrections to Bhabha scattering 
which is of crucial importance in determining the luminosity at the NLC (see the nice review of 
the current status in Rcf. [35]). Here there are 33 double box graphs to evaluate of which seven 
have been studied [36, 37, 38, 35]. Analytic expressions for the associated vertex graphs are also 
known [39, 35] and have been employed to calculate the QED [40] and QCD [41] corrections to the 
massive fermion form factor. 

In this paper, we take a first step towards calculating massless two-loop 2^2 scattering 
amplitudes with two off-shell legs. These processes include the NNLO QCD corrections to qq — ► 
V*V* (where V = W,Z) and the NLO corrections to gg — > V*V*. Here we indicate that the gauge 
bosons are off shell to account for resonance effects in the decay of the gauge boson. As in the on- 
shell and single off-shell cases, the MI include planar and non-planar box graphs as well as vertex 
and self energies. Altogether there are 11 planar box and 3 non-planar box master topologies, some 
of which require 2 or more MI to be computed. Our goal in this paper is much more modest. As 
a first step, we restrict ourselves to the planar and non-planar vertex graphs that form a simpler 
subset of the necessary MI. 

The MI presented here are ingredients for several interesting two-loop processes in their own 
right. In Higgs physics, the H — > V*V* decay receives QCD corrections when the Higgs couples 
to gluons (via a heavy top quark loop) which then couple to the weak bosons via a massless quark 
loop. This may be relevant for Higgs searches in the mass regions where the Higgs decays into 
two off-shell W bosons. In pure QCD, one can evaluate the two-loop triple gluon and quark-gluon 
two- loop vertices with massless quarks in a covariant gauge (as well as the gluon-ghost vertex). 1 
This is a useful input for Schwinger-Dyson studies of confinement as well as exploring how the 
Ward-Slavnov- Taylor identities generalise to the off-shell case. 

As in Refs. [15, 16, 39, 43, 44, 45, 35], we employ the differential equation technique to evaluate 
the MI. By building up from the simplest MI, we find differential equations in the external scales for 
new MI where the inhomogeneous terms involve known MI. Together with appropriate boundary 
conditions, the integral can then be systematically solved order by order in e. At each order we 
encounter one-dimensional integrals over the terms in the result for one order lower. These integrals 
yield polylogarithms and, because of the specific kinematics of the vertex graph with three off-shell 
legs, we find it necessary to extend the set of two-dimensional harmonic polylogarithms (2-d HPL) 
to include quadratic factors in the denominator (see also Ref. [45]). Our results are therefore 
presented in terms of the extended set of 2-d HPL's. 

Our paper is organised as follows. In Section 2, we give a brief summary of the differential 
equation method for solving MI while in Section 2. We define our notation and kinematics in 
Section 3 before we discuss harmonic polylogarithms (HPL) and the extended set of 2-d HPL 
necessary for these MI in Section 4. In any two-loop process, some of the graphs are factorisablc 

1 Note that Davydychev and Osland have studied the two-loop case where only one of the legs is off-shell [42] 
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and are given as products of one-loop graphs. Therefore, we list the one-loop MI in Section 5. 
We note that expressions to all orders in e for the one-loop triangle with off-shell legs have been 
provided in terms of log-sine functions by Davydychev [47] . The generalised log-sine functions can 
be directly related to Nielsen polylogarithms [48] and the all-order epsilon-expansion of one-loop 
massless vertex diagrams with three off-shell legs is given in terms of Nielsen polylogarithms in 
Ref. [48]. However, here we convert these results into 2-d HPL so that they can simply be combined 
with the other two-loop MI. Expressions for the two-loop MI are listed in Section 6, ordered by the 
number of propagators. Finally, our results are summarized in Section 7. 



2. Using differential equations to calculate master integrals 
2.1 The differential equations 

As discussed in Section 1, we use the differential equation (DE) method first suggested in [9] to 
evaluate the MPs. This method was expanded in detail in [46], and has since been used to calculate 
many two-loop MI, including those with multiple off-shell legs or internal masses [15, 16, 39, 43, 
44, 45, 35]. Here we present only a brief outline of this method, and for further details in the use 
of this method to find the e-expansion of master integrals we refer the reader to Refs. [15, 39]. 

The main idea of this method is to derive differential equations in external invariants for the 
master integrals. These equations are then solved using suitable boundary conditions to fix the 
constants of integration. 

The differential equations arc obtained by differentiating the MI with respect to the external 
momenta. Via the chain rule, linear combinations of these differentials are used to find the first 
order differential equations in the external invariants. These equations involve tensor loop integrals 
with additional powers of propagators. However, the IBP and LI identities can be used to simplify 
the equations 2 so that the differential equation is expressed only in terms of the MI itself and 
combinations of simpler topologies (i.e. integrals with less denominators). For this reason it is 
sensible to apply a 'bottom-up' approach, working from simpler to more complicated topologies, so 
that the MI is the only unknown in the differential equation. In other words, the inhomogeneous 
part of the differential equation is known (or at least the relevant terms in the e-expansion are 
known). For some topologies there is more than one MI leading to coupled differential equations. 
However as there is some freedom in the choice of which two-loop graphs to use as the master 
integrals, it simplifies the calculation to choose the master integrals such that they have different 
leading powers of e. In this case, the system of differential equations decouples on expansion in e 
(see, for example, [39]) . 

The differential equations are exact in the space-time dimension d, and can be solved as follows. 
Consider the inhomogeneous differential equation for the MI F with respect to the external scale 
x, 

±F(x)=A(x)F(x)+B(x). (2.1) 
If H(x) is a solution of the homogeneous equation 

^-H(x) = A(x)H(x). (2.2) 
ax 

then the full solution is given by 



F(x) = H(x) (^J H- 1 (x')B{x')dx' + Cj (2.3) 



2 To achieve this we have made extensive use of the Laporta algorithm [7]. We have coded our own version using 
FORM [49] and checked it against the Automatic Integral Reduction package AIR [50]. 



- 3 - 



where the constant C has to be fixed from the boundary conditions. These solutions are generally 
combinations of hypergeometric functions which arc difficult to expand in powers of e. Thus to 
find the e-expansion of the MI we must systematically expand each master integral F, and all 
rf-dependent terms of the differential equation in powers of e 



f = J2 A = E B = E 6V ( 2 - 4 ) 

i=—m i—0 i— — m 

where — m is the lowest power of e in the expansion and n is the highest power of e needed. It is 
assumed that the a 1 and b l are already known. Each coefficient /* satisfies the differential equation 
given by 

, m — i 

-f{x) = Y J ^{x)f- j {x)+b\x). (2.5) 

3=0 

It can be seen that the homogeneous part of all the equations generated by the e-expansion is simply 
the homogeneous solution H evaluated at d = 4 

h(x) = H(x)\ d=i (2.6) 

and so the solution is given by 

f(x) = h{x) ^ J" h-\x') «V).f ~V) + 6*(a/)^ dx' + c^j (2.7) 

where the constants c l have to be fixed from the boundary conditions at each order in e. Note 
that in general each coefficient /* will depend on p~ x so we solve the system of equations order 
by order, using repeated integration of the lower order results. It is for this reason that we require 
that A(x) has no poles in e as then /* would depend on f l+1 and the bottom-up approach would 
not be valid. 

2.2 The boundary conditions 

In general, the lowest order coefficient in e is determined solely by the boundary conditions. The 
boundary conditions are either obtained from the differential equation or from the master integral 
itself. To obtain limits from the differential equation it is necessary to examine the singular points 
in the coefficients of the differential equation. For example, if eq. 2.1 were to take the form 

JL F (x) - — !— F(x) + -J— B(x) (2.8) 

dx x — a x — a 

then we could multiply the whole equation by (a; — a) and let x — > a, then we have 

= F{x)\ x=a + B(x)\ x=a (2.9) 

giving the boundary condition on F(x). To obtain boundary conditions from the integral itself, we 
can use limits where the e expansion is known, for example where an offshell leg becomes massless. 
In both methods care has to be taken that the integral has a smooth limit at the chosen point so 
as not to miss or introduce hidden singularities. 

If H(x) is divergent at the boundary then the constant C is already determined by the necessary 
condition, 



H-\x')B(x')dx' + C = , (2.10) 
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which can be fulfilled by choosing the boundary point as the lower integration limit. The solution 
of the differential equation is then given by 



F(x) = H(x)F(x) = H(x) (^J H- 1 {x')B{x')dx^j . 



(2.11) 



It can be easily shown that this function satisfies the boundary condition. 
3. Kinematics 

We consider the vertex graph with three off-shell legs such that, 

-> Pi +P2+P3- 

The three kinematic scales are p\ . It is convenient to use the dimensionless variables 

_p\ _pI 
pi pi 

together with the determinant of the 2x2 Gram matrix, 

A = \(x, y) = \J(l - x - y) 2 - 4xy = \/(x - x )(x - X\) 

with 

xo = (l-Vv) 2 > X! = (l + ^y) 2 . 
In the following the dependence of A, xq and x\ on x and y will be implicitly understood. 



(3.1) 



(3.2) 



(3.3) 




Figure 1: The phase space for the vertex graph with three off-shell legs. The shaded region corresponds 
to the case where p| > p\, p\. The solid line marks the boundary where A = 0. 

In this paper we choose the kinematical configuration to suit the case of a heavy/offshell particle 
with momentum p% decaying into two lighter particles with momenta p\,p2- In this case the 
kinematically accessible region is given by the inequality, 

r ,2 



Pi + sJpl < JpI (3.4) 
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or after introduction of the dimensionless parameters x and y 

Vx~+^y<l. (3.5) 

The allowed parameter space is depicted by the shaded region in figure 1. Other regions are simply 
obtained by relabelling. 



4. Harmonic polylogarithms 

Solving the differential equations by repeated integration immediately suggests that the results 
be given in terms of Harmonic Polylogarithms (HPL's), whose properties are defined by repeated 
integration. The HPL's were first introduced in [2] as extensions of Nielsens polylogarithms, and 
later extended to 2-dimensions by [3]. In this section we briefly review the properties of one- 
dimensional HPL's (1-d HPL) before detailing the extension of the 2-d HPL's necessary to describe 
the vertex graph with three off-shell legs. 

4.1 1-d HPL's 

The weight, w, of a one-dimensional HPL, H( b;x), is the number of dimensions of the vector of 
parameters b . This vector, along with the argument x, fully describes the HPL. 

4.1.1 1-d HPL's with w = 1 

The weight-1 HPL's are defined as follows, 



HY„.y, _ / Iq f( a ,x')dx', a e {1,-1} 
[ ' } -\j*f(a,x')dx>, a = 



(4.1) 



where, 



/(i;*) = (4.2) 

f(0;x) = -, (4.3) 

x 

/(-i;-) -T^- (4-4) 

Note that H (0; x) is defined differently to avoid the logarithmic singularity at x = 0. Thus we have 

H(l;x) = -ln(l- x), (4.5) 

H(0;x) = lnx, (4.6) 

H{-l;x) = ln(l + a;), (4.7) 

and 

d 

—H(a;x) = f(a;x) with a e {+1,0,-1}. (4.8) 
ox 

4.1.2 1-d HPL's with w > 1 

The higher weight HPL's are recursively defined by, 

H{a,~b;x)= f dx'f{a;x')H(b;x), (4.9) 
Jo 

H(0,...,0;x) = -^-rln w x. (4.10) 
w< 
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Note that only the HPL's with weight vectors comprising only O's are defined differently, and are 
integrated between 1 and x to avoid logarithmic singularities. All others involve integration from 
to x. Under differentiation, the weight is reduced by unity, 

^-H (a, ~b; x) = /(a; x)H(b; x). (4.11) 
We also find the following relation useful: 

H n (0; x) = n!ff ((^0; x) (4.12) 

n zeros 

4.2 Extension of the 2-d HPL's 

The 2-d HPL's were introduced in [3] as the logical extension of the 1-d HPL's. The common 
extension is the linear basis 

f(a,x) = -^-, (4.13) 

f(-a,x) = ^±- x (4.14) 

However, while it is possible to use this basis to evaluate the integrals under investigation here, a 
more natural extension involves the square roots that are generated by Eq. 3.3. To this end, we 
extend the 2-d basis by the (quadratic) functions, 

f(X,x) = \, (4.15) 



A ' 

f(x\,x) = ±, (4.16) 



f(x ,x) = — , (4.17) 

x — xo 

f(x u x) = — . (4.18) 

x — x\ 

These functions are two dimensional, with explicit dependence on x and the dependence on y coming 
from x (y) and X\{y). The functions are chosen to be positive in the region, 

< x < 1, < y < 1, < A 2 , 

which immediately implies 

0< a; < rro = (1 - \/y) 2 < 1, < y < y = (1 - V^f < L ^+^y<l, (4.19) 
and which corresponds to the kinetically allowed region depicted in fig l. 3 



3 It is also possible to define a linear basis by making the Euler transformation, 

(t - a)(t + a) 

\ = 2t + x x = -- - 

(t-b) 

with 

1-2/ , 1 + 2/ 

a = , b = 

2 2 

such that 

dx dt dx dt 



A t-b A 2 (t - ai )(t - a 2 ) 

with 

=~(1 + VV) 2 , a 2 = -i(l- v ^) 2 - 
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Wc define the extended harmonic polylogarithms in the following way, 

G(a,w;x,y) = f f(a,x')G(w;x',y)dx', (4.20) 

J x 

and where the dependence on y is made explicit. Note that the lower boundary is x = x . The 
only exception from this definition is 

G(x ,...,x Q ;x,y)= f(xo,x')G(x ,...,x ;x',y)dx'. (4.21) 



However, HPL's of this form do not appear in the results presented in this paper. 

The choice of the integration limits is governed by the kinematic boundaries. To be able to 
evaluate the extended HPL's for x — ► xo, the condition lim x ^ Xo G(w;x,y) = has to be fulfilled. 
This ensures that 

fX 

f(x ,x')G(v;x',y)dx' (4.22) 



is finite in the limit x — > Xo. 

Note that HPL's of the form G(0, w; x, y) and G(x\,w; x,y) are divergent in the limit that 
x — > 0. This reflects the fact that taking the massless limit and making the e-expansion does not 
necessarily commute. In the cases where the x — > limit is smooth, we find that the HPL's appear 
in the combination, 

AG(0, w; x, y) = G(0, w; x, y) - (1 - y)G(x\, w; x, y) (4.23) 

which is finite as x — > 0. 

For the generalised HPL's of weight 1 we find, 

G{Q;x,y)=\og(^j, (4.24) 
G(A; x, y) = log ( ^^"^ - ( 4 - 25 ) 

G(x ;x,y) = -log(^— -^j , (4.27) 
G(x 1 -x 1 y) = -\og(^—^). (4.28) 

\Xi - X J 

In the course of solving the differential equations for the MI, we also find integrands of the 
form, 

i. 7 1 -—r^(x-x o r(x-x 1 f. 

\ z [x — Xo,i)A 

These are not independent and can be reduced to the set of quadratic basis functions using the 
following relation obtained via integration by parts, 



/ 



(x — x ) a (x — xi) /3 G(w, w; x, y)dx 

a + (3 + 2 
{P+l)(x - Xl ) 

+ , a , ttt 1 r / (x - x ) a+1 (x - Xi y +1 f{v)G{w- x, y)dx 



j (x — x ) a (x — xi) /3+1 G(w, w; x, y)dx 



{(3 + l)(x - x±) 
(x-x ) a+1 (x-x 1 f+ 1 

{{3 + 1)(X - X!) 



G(v,w,;x,y). (4.29) 
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5. One Loop 

All one-loop master integrals contain an overall factor 

(4^) e r(l + e)r 2 (l-e) 



5.1 Two point integrals 



16tt 2 T(l - 2e) 



(5.1) 




Figure 2: The one-loop Master Integral, BB(g 2 ) 

The only two point integral at one-loop level is the bubble graph shown in figure 2 which can 
be easily calculated using Feynman parameters and is given here for the sake of completeness, 



5.2 Three point integrals 



BB(q 2 ) = tS D (-q 2 )-<- 



1 



e(l-2e) 



(5.2) 



m% 




Figure 3: The one-loop Master Integral, Fo(mf,mJ,m3) 



The three point integrals with less than three massive legs are reducible to bubble integrals. 
The only master integral is the triangle with three massive legs, as shown in figure 3. The finite 
part of this diagram has been calculated in [51]. The (i-dimensional result for this integral can be 
found in [8, 52, 53], where it is given in terms of Appell functions. Davydychev calculated this 
integral to all orders in terms of log-sine functions [47] . 

In two-loop calculations this integral appears in products with other one-loop integral. To be 
able to combine these integrals with the genuine two loop integrals, it is necessary to express all 
integrals in terms of the same set of functions. To achieve this we apply the differential equation 
method also to the one loop triangle. 

The DE for Fo(m 2 , m^m 2 ) is given by 

m?A2 Hf = ((^r) A2 + (3 _ d)ml{ml ~ ml ~ ml) ) Fo 

+ (d - 3) (ml + m\- m^)BB(m 2 ) 
+ (d - 3) {ml + m 2 2 - ml) BB(to 2 ) 

- 2(d- 3)m 2 BB(m 2 ), (5.3) 

where, 

A 2 = m\ + m\ + m\ — 2m\m\ — 2m\m\ — 2m\m\. (5.4) 

The scalar triangle is completely symmetric under the interchange of the external scales. We choose 
to solve it in the configuration mf — pf with p\ > p 2 , p\. In this case, A 2 = p|A 2 (with A given 
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A 



j. To fix the constant 



by Eq. 3.3) and m\ — p 3 x, m 2 = p\y- By performing a transformation of variable m\ 
it can easily be seen that the homogeneous solution at d — 4 is Fq "' 
of integration it is necessary to look for suitable boundary points. Note that the limit x — > 
is not allowed because the one-loop with two external masses is divergent at d = 4. The only 
remaining possibility is to choose a point on the parabola given by A = 0, as shown in fig 1, i.e. 
x — ► £o,i = (1 ± y/y) 2 - I n this case the homogeneous solution is divergent at the boundary and so 
we apply the treatment discussed in section 2.3. The boundary condition for x — > x , corresponding 
to m\ — ► (m 2 — m 3 ) 2 is given by, 



F ( 



BB(m 2 2 



(m 2 - m 3 ) m 3 



+ 



BB((m 2 



m 3 ) 2 ) 



m 2 TO3 



+ 



BB(m 3 2 ) 
m 2 (m 2 - m 3 ) ' 



(5.5) 



Solving the differential equation order by order in e, we find the following expansion in e, 



^o(p 2 ,P 2 ,P 3 ) = 



where, 




*(-P§)- 




/5 



■n 2 ' n 2 



e 2 + C(e 3 ) 



(5.6) 



/ °(x, y) = - 2 G(A, 0; x,y)-2 G(A; a;, y) if (0; x ) 

+ G(A; a;, y) H(0; y) + (-1 + y) G(xA; a;, y) if (0; y), 



/o 1 ^, y) = + 2 G(0, A, 0; x, y) + 2 G(A, 0, 0; x, y) + 2 G(x , A, 0; x, y) 

+ 2 G( Xl , A, 0; x, y) + 2 G(A; x, y) H(0, 0; x ) - (G(A; x, y) H(0, 0; y)) 

+ (1 - y) G(xA; x, y) H(0, 0; y) + 2 G(0, A; x, y) H(0; x ) - (G(0, A; x, y) H(0; y)) 

+ (1 - y) G(0, xA; x, y) H(0; y) + 2 G(A, 0; x, y) H(0; x ) + 2 G(x , A; x, y) H(0; x ) 

- (G(x Q , A; x, y) H(0; y)) + (1 - y) G(x a ,x\; x, y) H(0; y) + 2 G(x u A; x, y) iJ(0; x ) 

- (G( Xl , A; x, y) H(0; y)) + (1 - y) G( Xl , xA; x, y) if (0; y), 

.f 2 (x, y) = - 2 G(0, 0, A, 0; x, y) - 2 G(0, A, 0, 0; x, y) - 2 G(0, x , A, 0; x, y) 

- 2 G(0, xi, A, 0; x, y) - 2 G(A, 0, 0, 0; x,y)-2 G(x , 0, A, 0; x, y) 

- 2 G(x , A, 0, 0; x, y) - 2 G(x , x , A, 0; x, y) - 2 G(x , x x , A, 0; x, y) 

- 2 G(xx , 0, A, 0; x, y) - 2 G(x x , A, 0, 0; x, y) - 2 G(xi,x , A, 0; x, y) 

- 2G(xi,xi, A,0;x,y) — 2 G(A; x, y) iJ(0, 0, 0; x ) + G(A; x, y) iJ(0, 0, 0; y) 

+ (-1 + y) G(xA; x, y) (0, 0, 0; y) - 2 G(0, A; x, y) ff(0, 0; x ) + G(0, A; x, y) if (0, 0; y) 
+ (-1 + y) G(0, xA; x, y) if(0, 0; y) - 2 G(A, 0; x, y) (0, 0; x ) - 2 G(x , A; x, y) if (0, 0; x ) 
+ G(x , A; x, y) if (0, 0; y) + (-1 + y) G(x , xA; x, y) if (0, 0; y) - 2 G(x 1; A; x, y) if (0, 0; x ) 
+ G(xi, A; x, y) if (0, 0; y) + (-1 + y) G(xi, xA; x, y) if (0, 0; y) - 2 G(0, 0, A; x, y) if (0; x ) 
+ G(0, 0, A; x, y) ff (0; y) + (-1 + y) G(0, 0, xA; x, y) if (0; y) - 2 G(0, A, 0; x, y) if (0; x ) 

- 2 G(0, x , A; x, y) ff (0; x ) + G(0, x , A; x, y) ff (0; y) + (-1 + y) G(0, x , xA; x, y) ff (0; y) 

- 2 G(0, xi, A; x, y) ff (0; x ) + G(0, x u A; x, y) ff (0; y) + (-1 + y) G(0, xi, xA; x, y) ff (0; y) 

- 2 G(A, 0, 0; x, y) ff (0; x ) - 2 G(x , 0, A; x, y) ff (0; x ) + G(x , 0, A; x, y) ff (0; y) 



- 10 - 



+ (-1 + y) G{x o ,0, x\- x, y) H(0; y) - 2 G(x , A, 0; x, y) H(0; x ) - 2 G(x , x , A; x, y) H(0; x ) 
+ G(x Ql x ,X;x,y) H(0;y) + (-1 + y) G(x , x , xX; x, y) H(0; y) - 2G(x ,x 1 ,X;x 1 y) H(0;x o ) 
+ G(x ai xi,X;x, y) H(0;y) + (-1 + y) G(xo, x\, xX; x,y) H(0; y) - 2 G(x\, 0, A;x, y) H(0:x o ) 
+ G{x u 0, A; x, y) H(0; y) + (-1 + y) G(x u O, xX; x, y) H(0; y) - 2 G(x u A, 0; x, y) H(0; x ) 
- 2 G(x!,x , A; x, y) H(0; x ) + G(xi, x , A; x, y) iJ(0; y) + (-1 + y) G(xi,xq, xX; x, y) H(0; y) 
-2G(x 1 ,x 1 , X;x,y) H(0;x o ) + G(x 1 ,xi 7 X;x,y) H(0:y) + (-1 + y) G(x 1 , x\, xX; x, y) H(0; y). 

We have checked that the e° term agrees with the results of Ref. [51] while the 0(e) term numerically 
agrees with that given in Ref. [47] . 



6. Two Loop 

6.1 Two point integrals 

There are only two two-point integrals, both of which can be obtained by repeated one-loop inte- 
gration and are widely available in the literature. We quote them here for the sake of completeness. 




Figure 4: The two- loop Master Integral, SS(<? 2 ). 
The three-propagator sunset graph SS(<? 2 ) is shown in Fig. 4 and is given by, 



r(3-3e)r(l-e)r(l + e) 



2\l-2e 



1 



13 115 
~8~ ~ To" e 



865 
"32" 



7 C(3) 



0(e 3 ) 



(6.1) 



'-OO- 

Figure 5: The two-loop Master Integral, GL(q 2 ). 
We denote the four-propagator graph shown in Fig. 6.1 by GL(q 2 ). It is given by, 

GL(<? 2 ) = BB(q 2 ) 2 = -i-fr 2 'Sl e2{l l _ 2e)2 

= -(-q 2 )- 2£ S 2 D ^ (1 + 4e + 12e 2 + 32e 3 + 80e 4 + G(e 5 )) . (6.2) 

6.2 Three point integrals 

6.2.1 Master Integrals with four propagators 

The simplest graph with four propagators is denoted by TGL and is shown in figure 6.2.1 and is 
simply the product of two bubbles, 

TGL(pIp 2 ) = BB(p 2 )BB(p 2 ) = -{-Plr^-PD-'SI l _ (6.3) 
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Figure 7: The two-loop Master Integrals, Fi(m\, m|, m§), ^(mi, ml, m§) 



There are two genuine two loop integrals with four propagators. These are denoted by F\ and 
F2 and shown in figure 6.2.1 The e expansion of these integrals is given to order e° in [54]. The 
corresponding two scale diagrams, when one of the legs becomes massless, are given in [15] to order 
e 2 in terms of HPL's. 

We obtain two coupled differential equations for F\ and F2, 

d 2 2 2 {d- A) m% (m\ -ml + ml) (d - 4) (-to? + m\ + to 2 ) 
f Fi(m 1I m 2 ,m 3 ) = — -, — rr-^ F 2 — F 1 



dm\ LK 17 27 A ' 2 (d-3) A 2 ' 2 A 2 

(3d -8) (-to 2 - to 2 + to 2 ) 2 , , (3d-8) go/ 2 , fRi , 
+ 2 A 2 to 2 88(771!) + — ^ — SS(m 2 ), (6.4) 



d F ,2 , m 2, - ((d- 3) (3 d - 10) (to 2 - to 2 + to 2 )) 
I<2{m 1 ,m2,m 3 ) 2 *i 



9to 2 zv 2 ' d/ 2 A 2 to 2 

(3d - 10)to 2 (to 2 - to 2 - to 2 ) - 2(d - 4)A 2 ) ) F 2 



2 A 2 to 2 

(d-3) (3d- 8) (3d-10) 2 ^ 
+ (d-4) A 2 to 2 bb(mi) 

(d-3) (3 d - 8) (3d - 10) (m 2 +m%- mf) 
2 (d - 4) A 2 to 2 to 2 

The homogeneous solutions at d = 4 are found to be 

j^hom -1 Tphom 
±> 1 = 1, ^2 



SS(to 2 ). (6.5) 



1 

A' 



Finding suitable boundary conditions for these coupled equations is more difficult than for the one 
loop triangle. Only F\ has a smooth limit for x — > (F 2 develops additional poles in e in that 
limit). On the other hand, taking the limit A — > in either differential equation only provides a 
single relation between the two integrals. Therefore we use the limit x — > to fix the integration 
constant of F\, matching the integral in this limit to the two-scale result given in [15]. Then we 
take the limit A — > (i.e. x — > xq) of F\ and use the limiting relation obtained from the differential 
equation to find F 2 in the limit A — > 0, thereby fixing the integration constant of F 2 . 
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The limiting relation in the limit x — > xq is, 

/ (d-3) ^i(wf,m|,m^)| m 2 =(m2 _ m3) 2 



J F 2 (TO 1 ,m 2 ,m 3 )| m 2 =(m2 _ m3)2 



m 2 m 3 



(d - 3) (3 d - 8) SS(m 2 2 ) (d - 3) (3 d - 8) SS(m§) 



(6.6) 



(d — 4) 77i 2 2 (m 2 — m 3 ) m 3 (d — 4) m 2 m 2 (m 2 — m 3 ) 



Without the second boundary constraint on F\ as x — > 0, this would only be sufficient to yield an 
independent determination of /{~ 2 only. 

The expansions in e are needed for two separate configurations. First, when the largest scale 
p\ is situated opposite the bubble (corresponding to m\ = p\) and second, when it lies adjacent to 
the bubble (m 2 = p\). For the first momentum configuration, we find that the expansions are as 
follows, 



where, 




sU-pIY 



E 

:-2,... 



2 2 

Pi P2 



fi\x,y) 
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-1 + x + y 
A 



G(A,0;x,y) + G(A,0;0, 2/ ) 



- \ AG(0; 0, y) H(0; y) + ^ G(0; x, y) H(0; y) 
+ -^ 1 + y) }: 1 + X + ^ G { xX;x,y)H { 0;y) 



+ 



2A 

G(X;0,y)H(0;x o 
-{-l + x + y) 



G(X;0,y)H(0;y) 



2A 



G(A; x, y) H(0; y) + 1 + * + V G(A; x, y) H(0; x ), 
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2 



H~ + 2 C3 + ~W AG(0; °' y) " AG(0 ' A ' 0; °' y) 



+ g G(0; x, y) + ^ G(A; 0, y) + - ( ^ + y) G( A: ,-. ,,) 



+ 
+ 



(it 2 {-1 + y) {-l + x + y)) 



12 A 

5 (-l + x + y) 



12 A 

G{xX;x,y) + 5G(X,0;0,y) 



G(A, 0; x, y) + 3( ^ + V) G(0, A, 0; x, y) 



- G(A, 0, 0; 0, y) + 2( 1 ± - + v) G(A, 0, 0; x, y) 



--G(A,0,A;0,y) 



A 

G(A, 0, xA; 0, y) + 1 G(A, A, 0; 0, y) 
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I G(A, A, 0; x, y) + -i±£ G(A, xA, 0; 0, y) - G(x , A, 0; 0, y) 



-1 + x + y 
A 



G(x , X,0;x,y) - G{x 1 , A, 0; 0, y) 



' ' ^ G(xi , A, 0; x, y) + 3 ( * + v) G(x\, X, 0; x, y) 



+ ^- if (1; y) - 5 if (1, 0; y) + 2 if (1, 0, 0; y) - if (1, 1, 0; y) 

- | AG(0; 0, y) if (0; y) + AG(0; 0, y) if (0, 0; y) 

- \ AG(0; 0, y) if (1, 0; y) - J AG(0, 0; 0, y) if (0; y) 

- AG(0, A; 0, y) if (0; x ) + ^ AG(0, A; 0, y) if (0; y) 

+ AG(0, xA; 0, y) if (0; y) - \ G(0; x, y) G(A, 0; 0, y) 

+ | G(0; x, y) if (0; y) - G(0; x, y) if (0, 0; y) + 1 G(0; x, y) if (1, 0; y) 
+ 5 G(A; 0, y) H(0; x ) - | G(A; 0, y) if (0; y) - 2 G(A; 0, y) fi(0, 0; x ) 
+ G(A; 0, y) if (0, 0; y) + 1 G(A; 0, y) if (1, 0; y) 

+ "^/"^^ G ( A ; y) G ( A > °; ' v) + 5 + x + g(a ; x , y ) ff(o ; x ) 

+ ' 5 i ~ 1 2 + x X + V) G(\; x, y) if (0; y) + 2 tl+£±V} G(A; x, y) if (0, 0; x ) 

+ ~ 1 + * + y G(A; x, y) if (0, 0; y) + ~ 1+ 2 * X + V G(X; x, y) if (1, 0; y) 

+ ( ' 1 + y)( 2 ' A 1 + X + y) G(xA;x,y)G(A,0;0,y) 

+ ' 5 tl±Atl±l±y} G(xA; x, y) if (0; y) 

+ { - 1 + V) { - 1 + X + y) G(xX;x,y)H(0Ay) 

+ ' (( ' 1 + ?;) 2 ( A 1 + X + ?/)) G(xA; x, y) if (1, 0; y) 

- 1 G(0, 0; x, y) if (0; y) + ~ 3 + V) G (0, A; x, y) if (0; x ) 

+ 3( ' 1 4 + A a: + y) G(0,A;s,g)ff(0;t/) 

3 (-1 + 2/) (-1 + x + y) 
H — G(0, xA; x, y) ii(0; y) 

- | G(A, 0; 0, y) if (0; x ) - \ G(A, 0; 0, y) if (0; y) 
+ ~ 2 + X + G(A, 0; x, y) if (0; x ) 

+ ~ ( ~ 1 4 + A a: + ?/) G(X, 0; x, y) if (0; y) 

+ i G(A, A; 0, y) if (0; x ) - \ G(A, A; 0, y) if (0; y) 

- | G(A, A; x, y) if (0; x ) + - A G(A, A; x, y) if (0; y) 
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+ G(X, xX; 0, y) H(0; x ) + 3 ( j + v) G(A, xX; 0, y) H(0; y) 

+ 3 { ~] + V) G(X, xX; x, y) H(0; y) - (G(xq, A; 0, y) H(0; x )) 



+ -G(x o ,X;0,y)H(0;y) 



-1 + x + y 



G(x ,X;x,y) H(0;x ) 



+ + G(x ,X; x, y) H(0; y) + -i±i G(x , xX; 0, y) H(0; y) 



2X 

j-l + y) (-1 + x + y) 
2X 



G(x ,xX; x, y) H(0; y) - (G( Xl , A; 0, y) H(0; x )) 



+ \ G{x\, A; 0, y) H(0; y) - ( 1 + * + y ) (7( A: .r. ,,)//((): ) 



+ l \l + V G(x 1 , A; x, y) H(0; y) + -i±£ 
, (-1 + 2/) {-1 + x + y) 

+ — G{xi,xX;x,y)H(0;y) 



G(x 1 ,xX;0,y) H(0;y) 



2X 

j-l+j/) (-1 + x + y) 

4A 
3 (-1 + 2/) 



GOrA,0;z,y)#(0;t/) 



G(xA, A;x,y)iJ(0;x o ) 



+ 



-3 (-1 + 2/) 



G(:cA,A;x,y)ff(0;j/) + 



-3(-l + 2/f 



G(xA,a;A;x,y)7J(0;7/) 



+ - A AG(0; 0, y) G(0; x, y) H(0; y) + 1 + V AG(0; 0, y) G(A; x, y) H(0; y) 

+ ' (( ' 1 + ?y) 4 ( A 1 + x + y)) AG (°; o. ») G ^ *. *) ^ *) 

- l - G(0; x, y) G(A; 0, y) H(0; x ) + \ G(0; x, y) G(A; 0, y) H(0; y) 



(-1 + x + y) 



+ 
+ 



2A 

-1 + x + y 



G(X;0,y)G(X;x,y)H(0;x o ) 



4A 

(-1 + y) (-1 + x + y) 
2X 



G(X;0,y)G(X;x,y)H(0;y) 



G(X;0,y)G(xX;x,y)H(0;x o ) 



and, 



*?(p?,J&Pl) = 



where, 




2\-l-2e 



Z^ 1 ^, y) = - 2 G(A, 0; x,y)-2 G(A; or, y) H(0; x ) 

+ G(A; I, y) H(0; y) + (-l + y) G(xX; x, y) H(0; y), 
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/ 2 °(x, y) = + -|- G(A; x, y) + - ( ^ + y) G^A; x, y) + 3 G(0, A, 0; x, y) 

+ 4 G(A, 0, 0; x, y) + 2 G(x , A, 0; x, y) + 2 G(x l7 A, 0; x, y) 

+ G(A; x, y) G(A, 0; 0, y) + 4 G(A; x, y) if (0, 0; x ) - 2 G(A; x, y) if (0, 0; y) 

- (G(A; x, y) if (1, 0; y)) + (1 - y) G(xA; x, y) G(A, 0; 0, y) 

+ (2 - 2 y) G(xA; x, y) #(0, 0; y) + (-1 + y) G(xA; x, y) #(1, 0; y) 

+ 3 G(0, A; x, y) if (0; x ) - | G(0, A; x, y) if (0; y) 

+ ~ 3 ^ + V) G(0, xA; x, y) if (0; y) + 4 G(A, 0; x, y) if (0; x ) 

+ \ G(A, 0; x, y) if (0; y) + 2 G(x , A; x, y) JT(0; x ) - (G(x , A; x, y) #(0; y)) 

+ (1 - y) G(x , xA; x, y) if (0; y) + 2 G(xi , A; x, y) if (0; x ) - (G(xi , A; x, y) #(0; y)) 

+ (1 - y) G(x! , xA; x, y) if (0; y) + G(xA, 0; x, y) if (0; y) 

- \ AG(0; 0, y) G(A; x, y) if (0; y) + AG ( ; 0, y) G(xA; x, y) if (0; y) 
+ G(A; 0, y) G(A; x, y) if (0; x ) - ± G(A; 0, y) G(A; x, y) if (0; y) 

+ (1 - y) G(A; 0, y) G(xA; x, y) if (0; x ) + G(A; 0, y) G(xA; x, y) if (0; y), 



7T 2 

/2 {x, y)=-( 3 G(A; x, y) + (-1 + y) Cs G(xA; x, y) + — G(0, A; x, y) 

+ + ^ G(0, xA; x, y) + ^ G(A, 0; x, y) + ^ G(x , A; x, y) 

-(tt 2 (-1 + y)) x X 

H g G(x ,xA;x,y) + y G(xi, A; x, y) 

-(tt 2 (-1 + y)) , -(^(-l + y)) 
+ — g ^G(xi,xA;x,y) + — ^ — ^G(xA,0;x,y) 

- I G(0, 0, A, 0; x, y) - 6 G(0, A, 0, 0; x, y) - 3 G(0, x , A, 0; x, y) 

- 3 G(0, xi, A, 0; x, y) - 8 G(A, 0, 0, 0; x, y) - | G(A, A, A, 0; x, y) 

+ 3 ^\ + V) G(X, xA, A, 0; x, y) - 3 G(x , 0, A, 0; x, y) - 4 G(x , A, 0, 0; x, y) 

- 2 G(x , x , A, 0; x, y) - 2 G(x , xi, A, 0; x, y) - 3 G(xi, 0, A, 0; x, y) 

- 4G(xi, A, 0,0; x,y) - 2 G(xi, x , A, 0; x, y) - 2G(x 1 ,x 1 , A,0;x,y) 

+ 3 { ~\ + V) G(x\ A, A, 0; x, y) + ~ 3 + ?j) G(xA, xA, A, 0; x, y) 

+ ^- AG(0; 0, y) G(A; x, y) + - { \\ + V) AG(0; 0, y) G(xA; x, y) 

- (AG(0, A, 0; 0, y) G(A; x, y)) + (-1 + y) AG(0, A, 0; 0, y) G(xA; x, y) 

+ ^ G(A; 0, y) G(A; x, y) + = ^ ^ + ^ G(A; 0, y) G(xA; x, y) 



(G(A; x, y) G(A, 0, 0; 0, y)) - \ G(A; x, y) G(A, 0, A; 0, y) 
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+ G(A; x, y) G(X, 0, xA; 0, y) + l - G(A; x, y) G(A, A, 0; 0, y) 

+ G(A; x, y) G(A, xX, 0; 0, y) - (G(A; x, y) G(x , A, 0; 0, y)) 

- (G(A; x, y) G{x u X, 0; 0, y)) + G(A; x, y) H(l; y) 

6 

- 8 G(A; x, y) if (0, 0, 0; x ) + 4 G(A; x, y) if (0, 0, 0; y) 
+ 2 G(A; x, y) if (1, 0, 0; y) - (G(A; x, y) if (1, 1, 0; y)) 

+ (-1 + y) G(xA; x, y) G(A, 0, 0; 0, y) + G(xA; x, y) G(A, 0, A; 0, y) 

+ ~ (-1 2 +y) G(xA; x, y) G(A, 0, xX; 0, y) + ^ G(xA; x, y) G(A, A, 0; 0, y) 
+ ~ (-1 2 +?/) G(xA; x, y) G(A, xA, 0; 0, y) + (-1 + y) G(xA; x, y) G(x , A, 0; 0, y) 

+ (-1 + y) G(xA; x, y) G(x 1; A, 0; 0, y) + ^ ^ + y) G(xA; x, y) if (1; y) 
+ 4 (-1 + y) G(xA; x, y) if (0, 0, 0; y) + (2 - 2 y) G(xA; x, y) JT(1, 0, 0; y) 
+ (-1 + y) G(xA; x, y) if (1, 1, 0; y) - | G(0, A; x, y) G(A, 0; 0, y) 

- 6 G(0, A; x, y) if (0, 0; x ) + 3 G(0, A; x, y) if (0, 0; y) 

+ | G(0, A; x, y) if (1, 0; y) + 3 ( ~ 2 +y) G(0, xA; x, y) G(A, 0; 0, y) 

+ 3 (-1 + y) G(0, xA; x, y) if (0, 0; y) + ~ 3 ^ + y) G(0, xA; x, y) if (1, 0; y) 

- 1 G(A, 0; 0, y) G(A, 0; x, y) - (G(A, 0; 0, y) G(x , A; x, y)) 

+ (-1 + y) G(A, 0; 0, y) G(x , xA; x, y) - (G(A, 0; 0, y) G(x 1; A; x, y)) 

— 1 + v 

+ (-1 + y) G(A, 0; 0, y) G(x l5 xX; x, y) + — ^ G(A, 0; 0, y) G(xA, 0; x, y) 

- 8 G(A, 0; x, y) if (0, 0; x ) - (G(A, 0; x, y) if (0, 0; y)) 

+ 1 G(A, 0; x, y) #(1, 0; y) - 4 G(x , A; x, y) if (0, 0; x ) 

+ 2 G(x , A; x, y) if (0, 0; y) + G(x , A; x, y) if (1, 0; y) 

+ 2 (-1 + y) G(x , xA; x, y) if (0, 0; y) + (1 - y) G(x , xA; x, y) if (1, 0; y) 

- 4 G(x! , A; x, y) if (0, 0; x ) + 2 G(xi , A; x, y) if (0, 0; y) + G(xi , A; x, y) if (1, 0; y) 
+ 2 (-1 + y) G(x! , xA; x, y) if (0, 0; y) + (1 - y) G(xj , xA; x, y) if (1, 0; y) 

+ (-1 + y) G(xA, 0; x, y) if (0, 0; y) + ±^ G(xA, 0; x, y) if (1, 0; y) 

- | G(0, 0, A; x, y) if (0; x ) + ^ G(0, 0, A; x, y) if (0; y) 
9 C— 1 4- wl 

+ V 4 ^ G(0, 0, xA; x, y) if (0; y) - 6 G(0, A, 0; x, y) if (0; x ) 

- G(0, A, 0; x, y) if (0; y) - 3 G(0, x , A; x, y) if (0; x ) 

+ | G(0, x , A; x, y) if (0; y) + 3 ( ~ 2 + y) G(0, x , xA; x, y) if (0; y) 

- 3 G(0, x l7 A; x, y) if (0; x ) + | G(0, an , A; x, y) if (0; y) 
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+ 3 ( ~2 + y) G(0, xi , xX; x, y) H(0; y) + 3 { ~\ + V) G(0, xX, 0; x, y) H(0; y) 

- 8 G(A, 0, 0; x, y) H(0; x ) - \ G(X, 0, 0; x, y) H(0; y) 

- | G(X, X, A; x, y) H(0; x ) + | G(A, A, A; x, y) ff(0; y) 

+ 3 ( ~^ + ^ G(X, X, xX; x, y) H(0; y) + 3 { ~\ + V) G(X, xX, A; x, y) H(0; x ) 

+ ~ 3 ( ~* + V) G(X, xX, A; x, y) H(0; y) + ~ 3 + y) G(A, xA, xA; x, y) if (0; y) 

- 3 G(x , 0, A; x, y) (0; x ) + | G(x , 0, A; x, y) (0; y) 

+ 3 + ^ G( - x °> °' xA; x ' y ) 4 G ( x °' A ' 0; x ' y ) x °) 

- ^ G(x , A, 0; x, y) iJ(0; y) - 2 G(x , x , A; x, y) H(0; x ) 

+ G(x , x , A; x, y) if(0; y) + (-1 + y) G(x , x , xA; x, y) ff(0; y) 

- 2 G(x , xi, A; x, y) i?(0; x ) + G(x , Xi, A; x, y) H (0; y) 

— 1 + y 

+ (-1 + y) G(x , xi, xA; x, y) iJ(0; y) H — G(x , xA, 0; x, y) iJ(0; y) 

- 3 G(xi , 0, A; x, y) if (0; x ) + | G(xi , 0, A; x, y) if (0; y) 

+ 3 + G(x! , 0, xA; x, y) if (0; y) - 4 G(*i , A, 0; x, y) if (0; x ) 

- ^G(xi, X,0;x,y)H(0;y) - 2 G(xi, x , A; x, y) H(0; x ) 

+ G(xi , x , A; x, y) if(0; y) + (-1 + y) G(xi , x , xA; x, y) ff(0; y) 

- 2 G(xi, xi, A; x, y) H(0;xo) 

+ G(xi , xi , A; x, y) if(0; y) + (-1 + y) G(xi , xi , xA; x, y) H(0; y) 

— 1 + v —1 + 77 

+ — G(x!, xA, 0; x, y) if (0; y) + — ^ G(xA, 0, 0; x, y) if (0; y) 

+ 3 + G ^ A > A ' A; 2/) # (°; *o) + ~ 3 + G(xA, A, A; x, y) H{0; y) 

+ ~ 3( ^ + j/) G(xA, A, xA; x, y) £T(0; y) + ~ 3 + V) G(xA, xA, A; x, y) ff(0; x ) 

+ 3 (~ 1 + 2/) G(xA, xA, A; x, y) H(0; y) + 3 (~ 1 + ^ G(xA, xA, xA; x, y) if (0; y) 

+ AG(0; 0, y) G(A; x, y) #(0, 0; y) - ^ AG(0; 0, y) G(A; x, y) Jf (1, 0; y) 

+ (1 - y) AG(0; 0, y) G(xA; x, y) Jf (0, 0; y) + AG ( ; 0, y) G(xA; x, y) H(l, 0; y) 

+ AG(0; 0, y) G(0, A; x, y) Jf (0; y) + ~ 3 ^ + V) AG(0; 0, y) G(0, xA; x, y) if (0; y) 

+ 1 AG(0; 0, y) G(A, 0; x, y) Jf (0; y) + l - AG(0; 0, y) G(x , A; x, y) if (0; y) 

+ ^ AG(0; 0, y) G(x , xA; x, y) if (0; y) + \ AG(0; 0, y) G(x 1; A; x, y) if (0; y) 

+ 1^ AG(0; 0, y) G(x l7 xA; x, y) Jf(0; y) + ^ AG(0; 0, y) G(xA, 0; x, y) Jf(0; y) 
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- 1 AG(0, 0; 0, y) G(A; x, y) 77 (0; y) + AG(0, 0; 0, y) G(xA; x, y) 77(0; y) 

- (AG(0, A; 0, y) G(A; x, y) 77(0; x )) + \ AG(0, A; 0, y) G(A; x, y) 77(0; y) 

+ (-l + y) AG(0, A; 0, y) G(xA; x, y) 77(0; x ) + ~ 3 ^ + y) AG(0, A; 0, y) G(xA; x, y) 77(0; y) 

+ AG(0, xA; 0, y) G(A; x, y) 77(0; y) + ztl+ll A G(0, xA; 0, y) G(xA; x, y) 77(0; y) 

- 2 G(A; 0, y) G(A; x, y) 77(0, 0; x ) + G(A; 0, y) G(A; x, y) 77(0, 0; y) 

+ X - G(A; 0, y) G(A; a;, y) 77(1, 0; y) + 2 (-1 + y) G(A; 0, y) G(xA; a:, y) 77(0, 0; x ) 
+ (1 - y) G(A; 0, y) G(xA; x, y) 77(0, 0; y) + ^ G(A; 0, y) G(xA; x, y) 77(1, 0; y) 

- | G(A; 0, y) G(0, A; x, y) 77(0; x ) + | G(A; 0, y) G(0, A; x, y) 77(0; y) 

+ 3 ^\ + V) G(A; 0, y) G(0, xA; x, y) 77(0; x ) + ~ 3 ^ + V) G(A; 0, y) G(0, xA; x, y) 77(0; y) 

- l - G(A; 0, y) G(A, 0; x, y) if (0; x ) + i G(A; 0, y) G(A, 0; x, y) 77(0; y) 

- (G(A; 0, y) G(x , A; x, y) #(0; x )) + l - G(A; 0, y) G(x , A; x, y) 77(0; y) 

+ (-1 + y) G(A; 0, y) G(x , xA; x, y) 77(0; x ) + ^ G(A; 0, y) G(x , xA; x, y) 77(0; y) 

- (G(A; 0, y) G( Xl , A; x, y) 77(0; x )) + ^ G(A; 0, y) G( Xl , A; x, y) if (0; y) 

+ (-1 + y) G(A; 0, y) G(x 1; xA; x, y) #(0; x ) + G(A; 0, y) G(xi, xA; x, y) 77(0; y) 

+ G(A; 0, y) G(xA, 0; x, y) 71(0; x ) + ^ G(A; 0, y) G(xA, 0; x, y) if (0; y) 

- | G(A; x, y) G(A, 0; 0, y) if (0; x ) - J G(A; x, y) G(A, 0; 0, y) 71(0; y) 
+ 1 G(A; x, y) G(A, A; 0, y) if (0; x ) - \ G(A; x, y) G(A, A; 0, y) if (0; y) 

+ G(A; x, y) G(A, xA; 0, y) 77(0; x ) + ~ 3 ^ + V) G(A; x, y) G(A, xA; 0, y) if (0; y) 

- (G(A; x, y) G(x , A; 0, y) 77(0; x )) + ± G(A; x, y) G(x , A; 0, y) 77(0; y) 

+ G(A; x, y) G(x , xA; 0, y) 77(0; y) - (G(A; x, y) G(xj , A; 0, y) 77(0; x )) 

+ l - G(A; x, y) G(x! , A; 0, y) 77(0; y) + G(A; x, y) G(xj , xA; 0, y) 77(0; y) 

+ 3 + C(xA; x, y) G(A, 0; 0, y) 77(0; x ) + G(xA; x, y) G(A, 0; 0, y) 77(0; y) 

+ ^ G(xA; x, y) G(A, A; 0, y) 77(0; x ) + G(xA; x, y) G(A, A; 0, y) 77(0; y) 

+ ~ (-1 2 +y) G(xA; x, y) G(A, xA; 0, y) 77(0; x ) + + ^ G (xA; x, y) G(A, xA; 0, y) 77(0; y) 

+ (-1 + y) G(xA; x, y) G(x , A; 0, y) 77(0; x ) + ^ G(xA; x, y) G(x , A; 0, y) 77(0; y) 

+ G(xA; x, y) G(x , xA; 0, y) 77(0; y) + (-1 + y) G(xA; x, y) Gfo , A; 0, y) 77(0; x ) 
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+ G(xA; x, y) G(xi , A; 0, y) H(0; y) + ( ^ + V) G(xX; x, y) G( Xl , xX; 0, y) H(0; y). 

The second momentum configuration is denned by, 



fI{pIpIpI) = 




= sU-Pi 



2\-2e 



-2,...,1 



2 2 
Pi P2 



= s£(-j>!) 



-l-2( 



2 2 

Pi P2 

Pi' Pi 



e l + 0(e 2 ) 



Expressions for f{ 2 in this momentum configuration are similarly lengthy to those for F^ 2 an d can 
be obtained in computer readable form from the authors. In each case, we have checked that the 
leading contribution agrees with the results of Ref. [54] . 

Note that solving the differential equation for the crossed triangle requires the functions Ff 2 (x, y) 
which arc symmetric in x and y and in addition the functions F 1 b 2 (x,y) and F 1 b 2 (y,x). The func- 
tions Fi 2 (y,x) are, of course, in principle known. However, exchanging x and y puts them in a 
form that is not suited for further integration over x. To get them into a suitable form we therefore 
recalculate them directly from the differential equations. 

6.2.2 Master Integrals with five propagators 




Figure 8: The two-loop Master Integral, TB(mf, m|, m%) 

There arc two master integrals with five propagators. The first is a product of one-loop integrals 
and the second is a genuine two-loop integral. The first, denoted by TB, is shown in figure 8 and 
is given by, 

TB(plplpl) = BB(p2)F (plpipl). (6.7) 

The e expansion for this integral is straightforwardly obtained from eq. 5.6 and we do not show it 
here. 



Figure 9: The two-loop Master Integral, Fz(m\ , m|, m%) 




The second irreducible two-loop diagram is denoted by F3 and is shown in figure 9. It is finite 
in 4-dimensions and the leading contribution has been calculated in Ref. [54]. 
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The differential equation for F 3 is given by, 

8 2 2 2 (d - 3) (-mf + mj + mj) 2 (3 d - 10) (d - 3) 2 2 2^ 
^F8(m 1) m 2l m 3 ) = * p ^F 3 _ 4 A2 F 1 (m 1 ,m 2 ,m 3 ) 

2 (3d- 10) (d- 3) „ . 2 2 2x - m 2 + m D ^ , 2 2 2x 

" (d-4 A2 Fl (m ' ' w? ' m 2 ) + A ~~tJ L F 2 (m? , m 2 , m 2 ) 



2 



( TO i - m|+m|) 2 2 2 , , 4(d-3) ^ riT , , ., 
p F 2 (m 3 ,m 1 ,m 2 ) + ^— — TGL(m 3 ,m 2 ). (6.8) 

The solution of the homogeneous equation at d — 4 is given by F 3 om = j and we therefore take 
the boundary condition at A = corresponding to x — > xq, so that 

F 3 (mi ,m 2 2 , ml)\ mlHm2 _ m3 y = - - _ _ Fl (m 2 2 ,m 3 2 , (m 2 - m 3 ) 2 ) 

(3d- 10) „ x2 2 2^ 

■Fi((m 2 -m 3 ) ,m 2 ,m 3 ) 



(d - 4) (m 2 - m 3 ) m 3 



F 2 (m 2 2 ,m 3 2 ,(m 2 -m 3 ) 2 ) (6.9) 



(d - 3) (m 2 - m 3 ) 
m 2 



(d - 3) m 3 



F 2 ((m 2 -m 3 ) ,m 2 , m 3 ) 




+ (d-4)Z'-i3)^ TGL((m2 - m3)2 ' W32) - 

Once again, there are two distinct kinematic configurations depending on the position of the large 
scale p 3 . For the first momentum configuration, the first two terms in the e expansion are given by, 



F 3 a (pipipi) =- — <^ I = su-Pir^U E /3(44V+o(e 2 ) 

\i=0....,l ^ 3 



where, 

/ 3 ° = - 6 Cs G(A; y) + 6 G(A, A, A, 0; x, y) 

+ ^- G(A; 0, y) G(A; x, y) + 2 G(A; x, y) G(A, 0, A; 0, y) 

- 2 G(A; x, y) G(A, A, 0; 0, y) + y G(A; x, y) (0; y) 

+ 2 G(A; x, y) if (0, 1, 0; y) - 2 G(A; x, y) H(l, 0, 0; y) 
+ G(A, 0; x, y) if (0, 0; y) + 6 G(A, A, A; x, y) H(0; x ) 

- 3 G(A, A, A; x, y) H(0; y) + (3 - 3 y) G(A, A, xA; x, y) ff(0; y) 

- AG(0; 0, y) G(A; x, y) if (0, 0; y) - AG(0, A; 0, y) G(A; x, y) if (0; y) 

- G(A; 0, y) G(A; x, y) #(0, 0; y) - 2 G(A; 0, y) G(A; x, y) (1, 0; y) 

- 2 G(A; x, y) G(A, A; 0, y) H(0; x ) + G(A; x, y) G(A, A; 0, y) H(0; y) 
+ 2 (-1 + y) G(A; x, y) G(A, xA; 0, y) £T(0; y), 



-3tt 4 

J3 = + 2Q G(A; x, y) - 3 (3 G(A, 0; x, y) + 6 (3 G(x , A; x, y) + 6 (3 G{x u A; x, y) 
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-(tt 2 + -(tt 2 (-1 + j/)) 
+ V 2 — G(A, A, xA; x, y) + V \ ^ G(A, xA, A; x, y) 

- 3 G(A, 0, A, A, 0; x, y) - 9 G(X, A, 0, A, 0; x, y) - 12 G(A, A, A, 0, 0; x, y) 

- 6 G(A, A, xo, A, 0; x, y) - 6 G(A, A, xi, X, 0; x, y) - 6 G(x , A, A, A, 0; x, y) 

- 6 G(x a , A, A, A, 0; x, y) + 3 Ca AG(0; 0, 2/) G(A; x, y) + y AG(0, A; 0, y) G(A; x, y) 

- 3 AG(0, A, 0, A; 0, y) G(A; a:, y) - 6 AG(0, A, A, 0; 0, y) G(A; x, y) + C 3 G(A; 0, y) G(A; x, y) 

2 2 2 

+ ^- G(A; 0, y) G(A, 0; x, y) + ^- G(A; 0, 2/) G(x , A; x, y) + ^- G(A; 0, y) G(x x , A; x, y) 
o 3 3 

+ ^- G(A; x, 2/) G(A, 0; 0, y) + - ^ + v) G(A; x, y) G(A, xA; 0, y) + y G(A; x, y) G(x , A; 0, y) 

+ y G(A; x, 2/) G(x! , A; 0, y) - 4 G(A; x, y) G(A, 0, 0, A; 0, y) - (G(A; x, y) G(A, 0, A, 0; 0, y)) 

- 2 G(A; x, y) G(A, A, 0, 0; 0, y) + (3 - 3 y) G(A; x, y) G(A, A, 0, xA; 0, y) 
+ (3 - 3 y) G(A; x, y) G(A, A, xA, 0; 0, y) + 4 G(A; x, y) G(A, x , A, 0; 0, y) 
+ 4 G(A; x, y) G(A, xj , A, 0; 0, y) + (3 - 3 y) G(A; x, y) G(A, xA, A, 0; 0, y) 

- 2 G(A; x, y) G(x , A, 0, A; 0, y) - 4 G(A; x, y) G(x , A, A, 0; 0, y) - 2 G(A; x, y) G(xi , A, 0, A; 0, y) 

- 4 G(A; x, y) G(xi , A, A, 0; 0, y) + 8 ( 3 G(A; x, y) £T (0; y) + 6 Ca G(A; x, y) H(l; y) 

2 2 2 

+ y^- G(A; x, y) if (0, 0; y) + y G(A; x, y) JT(0, 1; y) + y G(A; x, y) 0; y) 

+ 2 G(A; x, y) (0, 0, 0, 0; y) - 2 G(A; x, y) £T(0, 0, 1, 0; y) - 2 G(A; x, y) iJ(0, 1, 0, 0; y) 
+ 2 G(A; x, y) £T (0, 1, 1, 0; y) + 6 G(A; x, y) (1, 0, 0, 0; y) + 2 G(A; x, y) H(l, 0, 1, 0; y) 

- 2 G(A; x, y) H(l, 1, 0, 0; y) + 3 (-1 + y) G(A, 0; 0, y) G(A, A, xA; x, y) 

+ 3 (-1 + y) G(A, 0; 0, y) G(A, xA, A; x, y) - (G(A, 0; x, y) G(A, 0, A; 0, y)) 

+ G(A, 0; x, y) G(A, A, 0; 0, y) + ^- G(A, 0; x, y) if (0; y) - 3 G(A, 0; x, y) if (0, 0, 0; y) 

- (G(A, 0; x, y) ff(0, 1, 0; y)) + G(A, 0; x, y) if (1, 0, 0; y) - 2 G(x , A; x, y) G(A, 0, A; 0, y) 

-TT 2 

+ 2 G(x , A; x, y) G(A, A, 0; 0, y) + — G(x , A; x, y) if (0; y) 

- 2 G(x , A; x, y) if (0, 1, 0; y) + 2 G(x , A; x, y) if (1, 0, 0; y) - 2 G(x 1; A; x, y) G(A, 0, A; 0, y) 

-TT 2 

+ 2 G(xi, A; x, y) G(A, A, 0; 0, y) + — G(x 1; A; x, y) if (0; y) - 2 G(x 1; A; x, y) #(0, 1, 0; y) 
+ 2 G(X! , A; x, y) if (1, 0, 0; y) - l - G(A, 0, 0; x, y) if (0, 0; y) - 12 G(A, A, A; x, y) #(0, 0; x ) 
+ H G(A, A, A; x, y) if (0, 0; y) + 9 + G(A, A, xA; x, y) if (0, 0; y) 
+ (3 - 3 y) G(A, A, xA; x, y) #(1, 0; y) + ~ 3 ^ + V) G(A, xA, A; x, y) if (0, 0; y) 

+ (3 - 3 y) G(A, xA, A; x, y) if (1, 0; y) + ~ 3 + y) G(A, xA, xA; x, y) (0, 0; y) 

- (G(xo, A, 0; x, y) if (0, 0; y)) - (G(xi, A, 0; x, y) if (0, 0; y)) - 3 G(A, 0, A, A; x, y) if (0; x ) 
+ | G(A, 0, A, A; x, y) /f (0; y) + 3 + G(A, 0, A, xA; x, y) if (0; y) 

- 9 G(A, A, 0, A; x, y) if (0; x ) + | G(A, A, 0, A; x, y) /f (0; y) 
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9 (— 1 + v) 

+ 1 2 ~ G(X, A, 0, xA; x, y) if (0; y) - 12 G(A, A, A, 0; x, y) ff (0; x ) 

- 6 G(A, A, x , A; a;, y) if (0; x ) + 3 G(A, A, x , A; a;, y) ff (0; y) 

+ 3 (-1 + y) G(A, A,x ,xA;x,y) H(0;y) — 6 G(A, A, x\, A; x, y) ff(0;x ) 
+ 3G(A,A,zi,A;a:,y)ff(0;y) + 3 (-1 + y) G(A, A, x l7 xA; x, y) ff(0; y) 

+ 3 G ( A ' A ' xA ' 0; x < f ) y ) + 3 ^2 + y) G(A ' XX > A ' ° ; X ' y) H ( ° ; y) 

- 6 G(x , A, A, A; x, y) ff(0; x ) + 3 G(x , A, A, A; x, y) if (0; y) 

+ 3 (-1 + y) G(x Q , A, A, xA; x, y) #(0; y) - 6 G(x l7 A, A, A; x, y) H(0; x ) 
+ 3G(x 1 ,A,A,A;x,y)ii(0;y) + 3 (-1 + y) G(xi, A, A, xA; x, y) H(0; y) 

+ ^- AG(0; 0, y) G(A; x, y) ff (0; y) + 3 AG(0; 0, y) G(A; x, y) ff (0, 0, 0; y) 
+ AG(0; 0, y) G(A; x, y) if (0, 1, 0; y) - (AG(0; 0, y) G(A; x, y) if (1, 0, 0; y)) 
+ i AG(0; 0, y) G(A, 0; x, y) ff (0, 0; y) + AG(0; 0, y) G(x , A; x, y) if (0, 0; y) 

+ AG(0; 0, y) G(x 1; A; x, y) if (0, 0; y) + ~ 3 ^ + V) AG(0; 0, y) G(A, A, xA; x, y) ff(0; y) 

+ ~ 3 ^ + V) AG(0; 0, y) G(A, xA, A; x, y) if(0; y) - 1 AG(0, 0; 0, y) G(A; x, y) ff (0, 0; y) 

- (AG(0, A; 0, y) G(A; x, y) G(A, 0; 0, y)) + | AG(0, A; 0, y) G(A; x, y) ff(0, 0; y) 

+ AG(0, A; 0, y) G(A; x, y) ff(l, 0; y) + 1 AG(0, A; 0, y) G(A, 0; x, y) ff(0; y) 
+ AG(0, A; 0, y) G(x , A; x, y) if (0; y) + AG(0, A; 0, y) G(x x , A; x, y) if (0; y) 
+ AG(0, xA; 0, y) G(A; x, y) fi(0, 0; y) + | AG(0, 0, A; 0, y) G(A; x, y) if (0; y) 

- i AG(0, A, 0; 0, y) G(A; x, y) ff (0; y) - 8 AG(0, A, A; 0, y) G(A; x, y) ff (0; x ) 

+ 4 AG(0, A, A; 0, y) G(A; x, y) ff (0; y) + 5 { ~\ + V) AG(0, A, xA; 0, y) G(A; x, y) ff (0; y) 
+ AG(0, x , A; 0, y) G(A; x, y) if (0; y) + AG(0, x l7 A; 0, y) G(A; x, y) ff(0; y) 

+ 1 g> AG(0, xA, A; 0, y) G(A; x, y) ff (0; y) + — G(A; 0, y) G(A; x, y) if(0; y) 
2 D 

+ ^- G(A; 0, y) G(A; x, y) ff (1; y) + 3 G(A; 0, y) G(A; x, y) ff (0, 0, 0; y) 

+ G(A; 0, y) G(A; x, y) ff (0, 1, 0; y) + 3 G(A; 0, y) G(A; x, y) ff(l, 0, 0; y) 

- 2 G(A; 0, y) G(A; x, y) (1, 1, 0; y) + \ G(A; 0, y) G(A, 0; x, y) ff(0, 0; y) 

+ G(A; 0, y) G(A, 0; x, y) ff(l, 0; y) + G(A; 0, y) G(x , A; x, y) ff(0, 0; y) 

+ 2 G(A; 0, y) G(x , A; x, y) ff(l, 0; y) + G(A; 0, y) G(x l7 A; x, y) ff(0, 0; y) 

+ 2 G(A; 0, y) G(xi, A; x, y) ff(l, 0; y) + 3 (-1 + y) G(A; 0, y) G(A, A, xA; x, y) ff(0; x ) 

+ ~ 3 ^ + y) G(A; 0, y) G(A, A, xA; x, y) ff(0; y) + 3 (-1 + y) G(A; 0, y) G(A, xA, A; x, y) ff (0; x ) 

+ ~ 3 ^ + V) G(A; 0, y) G(A, xA, A; x, y) ff (0; y) + G(A; x, y) G(A, 0; 0, y) 2 

+ (1 - y) G(A; x, y) G(A, 0; 0, y) G(A, xA; 0, y) - 2 G(A; x, y) G(A, 0; 0, y) G(x , A; 0, y) 
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- 2 G(A; x, y) G(X, 0; 0, y) G(x u X; Q,y)-\ G(A; x, y) G(A, 0; 0, y) H(0, 0; y) 

- (G(A; x, y) G(A, 0; 0, y) if (1, 0; y)) + 4 G(A; x, y) G(A, A; 0, y) H(0, 0; x ) 

- | G(A; x, y) G(A, A; 0, y) if (0, 0; y) + ~ ? ^ + V) G(A; x, y) G(A, xA; 0, y) if (0, 0; y) 

+ (-1 + y) G(A; x, y) G(A, xA; 0, y) Ji(l, 0; y) + G(A; x, y) G(x , A; 0, y) H(0, 0; y) 
+ 2 G(A; x, y) G(x , A; 0, y) 0; y) + (-1 + y) G(A; x, y) G(x , xA; 0, y) Jf (0, 0; y) 
+ G(A; x, y) G(an , A; 0, y) H(0, 0; y) + 2 G(A; x, y) G(an , A; 0, y) (1, 0; y) 
+ (-1 + y) G(A; x, y) G(x u xX; 0, y) JT(0, 0; y) - (G(A; x, y) G(A, 0, A; 0, y) Ji(0; x )) 

- | G(A; x, y) G(A, 0, A; 0, y) if (0; y) + (1 - y) G(A; x, y) G(A, 0, xA; 0, y) Jf (0; y) 
+ 3 G(A; x, y) G(A, A, 0; 0, y) Jf (0; x ) + i G(A; x, y) G(A, A, 0; 0, y) JJ(0; y) 

+ (5 - 5 y) G(A; x, y) G(A, A, xA; 0, y) if (0; x ) + 5 ^~\ +V) G(X; x, y) G(A, A, xA; 0, y) Jf (0; y) 

+ 2 G(A; x, y) G(A, x , A; 0, y) JJ(0; x ) - (G(A; x, y) G(A, x , A; 0, y) ff(0; y)) 

+ (2 - 2 y) G(A; x, y) G(A, x , xA; 0, y) Jf(0; y) + 2 G(A; x, y) G(A, Xi , A; 0, y) if (0; x ) 

- (G(A; x, y) G(A, x l7 A; 0, y) ff(0; y)) + (2 - 2 y) G(A; x, y) G(A, x l7 xA; 0, y) JJ(0; y) 

+ (-1 + y) G(A; x, y) G(A, xA, 0; 0, y) Jf (0; y) + (5 - 5 y) G(A; x, y) G(A, xA, A; 0, y) JJ(0; x ) 

+ 5 G ( A; f ) G ( A ' xX ' A; °' f ) H (° ; f ) + 4 (-1 + y) 2 G(A; x, y) G(A, xA, xA; 0, y) JJ(0; y) 

+ G(A; x, y) G(x , 0, A; 0, y) Jf (0; y) - 8 G(A; x, y) G(x , A, A; 0, y) if (0; x ) 

+ 4 G(A; x, y) G(x , A, A; 0, y) if (0; y) + 2 (-1 + y) G(A; x, y) G(x , A, xA; 0, y) JJ(0; y) 

+ 2 (-1 + y) G(A; x, y) G(x , xA, A; 0, y) #(0; y) + G(A; x, y) G(xj , 0, A; 0, y) Jf(0; y) 

- 8 G(A; x, y) G(x 1; A, A; 0, y) Jf(0; x ) + 4G(A; x, y) G(x 1; A, A; 0, y) JJ(0; y) 

+ 2 (-1 + y) G(X;x,y)G(x 1 ,X,xX;0,y) H(0;y) + 2 {-1 + y) G(X; x, y) G(x\ , xX, A; 0, y) Jf (0; y) 
+ G(A, 0; x, y) G(A, A; 0, y) Jf(0; x ) - \ G(A, 0; x, y) G(A, A; 0, y) JJ(0; y) 

+ (1 - y) G(A, 0; x, y) G(A, xA; 0, y) JJ(0; y) + 2 G(A, A; 0, y) G(x , A; x, y) #(0; x ) 

- (G(A, A; 0, y) G(x , A; x, y) Jf(0; y)) + 2 G(A, A; 0, y) G(xi, A; x, y) JJ(0; x ) 

- (G(A, A; 0, y) G(xj , A; x, y) Jf(0; y)) + (2 - 2 y) G(A, xA; 0, y) G(x , A; x, y) #(0; y) 
+ (2 - 2 y) G(A, xA; 0, y) G(xi , A; x, y) Jf(0; y). 

The other momentum configuration required is defined by, 




Expressions for /| in this momentum configuration arc similarly lengthy to those for Fg and can 
be obtained in computer readable form from the authors. In each case, we have checked that the 
finite contribution agrees with the results of Ref. [54] . 

6.2.3 Master Integrals with six propagators 

The only MI with 6 propagators is the crossed triangle which is denoted by F4 and is shown in 
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Figure 10: The two-loop Master Integral, i<4(mi, ra 2 , ml) 

figure 10. 

This integral is fully symmetric in all three legs, and satisfies the differential equation, 
d _ (6 - d) (- mi 2 + m 2 2 + m 3 2 ) 8 _ . 2 2 2 , 



dm 2 4 A 2 A 2 

4(m 1 2 +m 2 2 -m 3 2 ) 4 (-to,. 2 + m 2 2 - m 3 2 ) 

A 2 ^ 2 2(m2 ' m3 ' mi ) A 2 ^ 2 2(m3 ' m2 ' mi ^ 

4(d~4)i 3 (m 1 2 ,m 2 2 ,m 3 2 ) 2 (d - 4) ( mi 2 + m 2 2 - m 3 2 ) 2 2 

7^ 1 ~To 5 -f 3 m 2 ,m 3 >i ) 

A 2 A 2 mi 2 

2(d-4) (-m 1 2 +m 2 2 -m 3 2 ) 2 2 2 
A2mi2 L F 3 (m 3 ,m 2 ,mi ). (6.10) 

The homogeneous solution at = 4 is F^ 10 ™ = A -2 while the boundary condition at A = 
corresponding to x — > x is given by, 

F 4 (to 2 , to 2 , m 3 )| m 2 = ( m2 _ m3 )2 = y_ m 6 2 - m ^ 3 F 2 (rn|,rn 2 , (m 2 - ™ 3 ) 2 ) 



4(m 2 -m 3 ) 2 , 
(a — 6) to 2 m 3 (rf — 6) m\ m 2 



-F 2 ((m 2 - m 3 ) ,m 2 ,m 3 )--r- — — — 2 ^(m 2 , m 2 , (m 2 - m 3 ) ) 



i 2(d-4) (m 2 -m 3 ) 2 2 2 2(d-4) F m ^2 2 

H 71 — — 2 ^3\m 2 ,m 3 ,(m 2 - m 3 ) ) - — — — r 3 ((m 2 - m 3 J , m 2 ,m 3 j 

(d — 6) mj m 3 (a - 6) m 2 m 3 

2 (d - 4) (m 2 - m 3 ) 



-F 3 (m 3 ,m 2 , (m 2 - m 3 ) ). (6.11) 



(<i — 6) to 2 m 2 

We find that the first two terms of the e-expansion are given by, 




Pi ^7 I ^ /« 2 n 2 \ 



E fi( 

i=0,...,l v 



where, 



.0/ x -87T 2 , . -47T 2 (-1 + y) , . -47T 2 (-1 + w) x s 

/ 4 °(x, y) = + — ^— G(A, A; x, y) + K - ^ G(A, xA; x, y) + K - ^ G(xX, A; x, y) 

- 8 G(A, 0, A, 0; x, y) - 16 G(A, A, 0, 0; x, y) + 16 G(A, 0; 0, x ) G(A, A; x, y) 

+ 8 (-1 + y) G(A, 0; 0, y) G(A, xA; x, y) + 8 (-1 + y) G(A, 0; 0, y) G(xA, A; x, y) 

- 16 G(A, A; x, y) if (0, 0; x ) - 4 G(A, A; x, y) if (0, 0; y) - 16 G(A, A; x, y) if (1, 0; x ) 
+ (4 - 4 y) G(A, xA; x, y) if (0, 0; y) + (8 - 8 y) G(A, xA; x, y) if (1, 0; y) 

+ (4 - 4 y) G(xA, A; x, y) if (0, 0; y) + (8 - 8 y) G(xA, A; x, y) #(1, 0; y) 

- 4 (-1 + y) 2 G(xA, xA; x, y) if (0, 0; y) - 8 G(A, 0, A; x, y) H(0; x ) 
+ 4 G(A, 0, A; x, y) if (0; y) + 4 (-1 + y) G(A, 0, xA; x, y) if (0; y) 

- 16 G(A, A, 0; x, y) if (0; x ) + 8 G(A, A, 0; x, y) if (0; y) 
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+ 4 (-1 + y) G(A, xX, 0; x, y) H(0; y) + 4 (-1 + y) G(xX, X, 0; x, y) H{0; y) 

- 8 AG(0; 0, x ) G(A, A; x, y) H(0; x ) + (4 - 4 y) AG(0; 0, y) G(A, xA; x, y) iJ(0; y) 
+ (4 - 4 y) AG(0; 0, y) G(xA, A; x, y) £T (0; y) - 8 G(A; 0, x ) G(A, A; x, y) #(0; x ) 
+ 16 G(A; 0, x ) G(A, A; x, y) #(0; y) + 8 (-1 + y) G(A; 0, y) G(A, xA; x, y) £T(0; x ) 

+ (4 - 4y) G(A; 0, y) G(A, xA; x, y) #(0; y) + 8 (-1 + y) G(A; 0, y) G(xA, A; x, y) ff(0; x ) 
+ (4 - 4y) G(A; 0, y) G(xA, A; x, y) if (0; y) 
/i(x,y) = + 112C 3 G(A,A;x,y)+8 (-1 + y) C 3 G(X, xA; x, y) + 8 (-1 + y) C3 G(xA, A; x, y) 

+ 32 (-1 + y) 2 C.3 G(xA, xA; x,y) + ^- G(A, 0, A; x, y) + 4 ^ + G(A, 0, xA; x, y) 

4tt 2 ^ x „ \ 871- 2 , , 47r 2 (-1 + y) _ A . 
+ — G(A, A, 0; x, y) + -y- G(A, x , A; x, y) H G(A, x , xA; x, y) 

8tt 2 ,,,, x x 47r 2 (-l + y)^., , , 47r 2 (-l + y)^., 
+ — G(A, X! , A; x, y) + K — v -± G(A, x x , xA; x, y) + K — J -± G(A, xA, 0; x, y) 

-8 71-2 ^/ x x x -4tt 2 (-1 + y) ^. , -4tt 2 (-1 + y) x 
H — G(x , A, A; x, y) H G(x , A, xA; x, y) H G(x , xA, A; x, y) 

-8 71-2 ^/ xx x - 4 71-2 (-1 +2/) ^/ x x x -4tt 2 (-1 +y) , , 

H — G(xi,A, A;x,y) H G(xi, A, xA; x, y) H G(xi, xX, A; x, y) 

+ ( 3 1 + j/) G(xA, 0, A; x, y) + ^ + G(xA, A, 0; x, y) + ^ + G(xX, x , A; x, y) 

47T 2 f-l + w) -47T 2 (-1 + v) 2 

+ \ y! G(xA, Xi , A; x, y) + V ^ G(xA, xA, 0; x, y) + 16 G(A, 0, 0, A, 0; x, y) 

+ 24 G(A, 0, A, 0, 0; x, y) + 8 G(A, 0, x , A, 0; x, y) + 8 G(A, 0, x u \ 0; x, y) 

+ 48 G(A, A, 0, 0, 0; x, y) - 72 G(A, A, A, A, 0; x, y) + 8 G(A, x , 0, A, 0; x, y) 

+ 16 G(A, x , A, 0, 0; x, y) + 8 G(A, x u 0, A, 0; x, y) + 16 G(A, x u X, 0, 0; x, y) 

+ 24 (-1 + y) 2 G(A, xA, xA, A, 0; x, y) - 8 G(x , A, 0, A, 0; x, y) - 16 G(x , A, A, 0, 0; x, y) 

- 8G(xi, A,0, A,0;x,y) - 16 G(xi, A, A, 0, 0; x, y) + 24 (— 1 + y) 2 G(xA, A, xA, A, 0; x, y) 

— 2 7T 2 

+ 24 (-1 + y) 2 G(xA, xA, A, A, 0; x, y) + — ^— AG(0; 0, x ) G(A, A; x, y) 

+ 2?r2 ( ~ 3 1 + y) A G(0; 0, x ) G(xA, xA; x, y) + AG (°; °. 2/) G ( A > x > 2/) 

+ g'^^ AG (0; 0, y) G(xA, A; x, y) - 12 AG(0, A, 0; 0, x ) G(A, A; x, y) 

+ 12 (-1 + y) 2 AG(0, A, 0; 0, x ) G(xA, xA; x, y) - 16 (-1 + y) AG(0, A, 0; 0, y) G(A, xA; x, y) 

- 16 (-1 + y) AG(0, A, 0; 0, y) G(xA, A; x, y) + G(A; 0, x ) G(A, A; x, y) 

+ 27r2( ~ 1 + y) G(A; 0, x ) G(xA, xA; x, y) + ^ G(A; 0, y) G(A, A; x, y) 

- 8 G(A, 0; 0, x ) G(A, 0, A; x, y) - 8 G(A, 0; 0, x ) G(A, A, 0; x, y) 

- 16G(A,0;0,x )G(A,x , A;x,y) - 16 G(A, 0; 0, xo) G(A, x\ , A; x, y) 
+ 16G(A,0;0,x )G(x ,A, A;x,y) + 16 G(A, 0; 0, x ) G{x u \ A; x, y) 

+ 8 (-1 + yf G(A, 0; 0, x ) G(xA, xA, 0; x, y) + (8 - 8 y) G(A, 0; 0, y) G(A, 0, xA; x, y) 
+ (8 - 8 y) G(A, 0; 0, y) G(A, x , xA; x, y) + (8 - 8 y) G(A, 0; 0, y) G(A, x 1; xA; x, y) 
+ (8 - 8 y) G(A, 0; 0, y) G(A, xA, 0; x, y) + 8 (-1 + y) G(A, 0; 0, y) G(x , A, xA; x, y) 
+ 8 (-1 + y) G(A,0;0,y)G(x ,xA,A;x,y)+8 (-1 + y) G(A, 0; 0, y) G(x u A, xA; x, y) 
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+ 8 (-1 + y) G(A, 0; 0, y) G(xi, xA, A; x, y) - 16 (-1 + y) G(A, 0; 0, y) G(x\, 0, A; x, y) 
+ (8 - 8 y) G(A, 0; 0, y) G(xA, A, 0; x,y) + (8-8 y) G(X, 0; 0, y) G{xX, x , A; x, y) 
+ (8 - 8 y) G(X, 0; 0, y) G(xX, x x , A; x, y) - 16 G(X, A; x, y) G(X, 0, 0; 0, x ) 

- 16 G(X, A; x, y) G(A, 0, A; 0, y) - 12 G(A, A; x, y) G(A, A, 0; 0, x ) 
+ 16 G(A, A; x, y) G(A, A, 0; 0, y) - 8 G(A, A; x, y) G(x , A, 0; 0, x ) 

- 8 G(A, A; x, y) G( Xl , A, 0; 0, x ) + ^- G(A, A; x, y) if (0; x ) 

+ G(A, A; x, y) if (0; y) + ^f- G(A, A; x, y) if (1; x ) 

+ 48 G(A, A; x, y) if (0, 0, 0; x ) + 12 G(A, A; x, y) if (0, 0, 0; y) 
+ 8 G(A, A; x, y) if (0, 1, 0; x ) - 16 G(A, A; x, y) if (0, 1,0; y) 
+ 8 G(A, A; x, y) if (1, 0, 0; x ) + 16 G(A, A; x, y) if (1, 0, 0; y) 

- 8 G(A, A; x, y) if (1, 1, 0; ar ) - 16 (-1 + y) G(A, xX; x, y) G(A, 0, 0; 0, y) 

+ 8 (-1 + yf G(A, xX; x, y) G(A, 0, xX; 0, y) + 8 (-1 + y) 2 G(A, xX; x, y) G(A, xX, 0; 0, y) 
-16 (-l + y) G(X,xX;x,y)G(x ,X,0;0,y)-16 (-1 + y) G(X, xX; x, y) G(xi , X, 0; 0, y) 

+ 4?r2 ( 3 1 + ?/) G(A, xX; x, y) H(0: y) + t± ± y) G(A, xX; x, y) H(l;y) 

+ 12 (-1 + y) G(A, xA; x, y) H(0, 0, 0; y) + 8 (-1 + y) G(A, xA; x, y) if(0, 1, 0; y) 
+ 24 (-1 + y) G(A, xA; x, y) 0, 0; y) - 16 (-1 + y) G(A, xA; x, y) 1, 0; y) 

- 16 (-1 + y) G(xA, A; x, y) G(A, 0, 0; 0, y) + 8 (-1 + y) 2 G(xA, A; x, y) G(A, 0, xA; 0, y) 
+ 8 (-1 + y) 2 G(xA, A; x, y) G(A, xA, 0; 0, y) - 16 (-1 + y) G(xA, A; x, y) G(x , A, 0; 0, y) 

47T 2 f-1 + «1 

- 16 (-1 + y) G(xA, A; x, y) G(x x , A, 0; 0, y) + y - ^ G(xA, A; x, y) if (0; y) 

+ ~ 87f2 ( 3 ~ 1 + ?y) G (' xA ' y) y) + 12 (-1 + y) G(xA, A; x, y) if (0, 0, 0; y) 
+ 8 (-1 + y) G(xA, A; x, y) H(0, 1, 0; y) + 24 (-1 + y) G(xA, A; x, y) 0, 0; y) 

- 16 (-1 + y) G(xA, A; x, y) 1, 0; y) + 16 (-1 + yf G(xA, xA; x, y) G(A, 0, 0; 0, x ) 

+ 12 (-1 + yf G(xA, xA; x, y) G(A, A, 0; 0, x ) + 8 (-1 + yf G(xA, xA; x, y) G(x , A, 0; 0, x ) 

+ 8 (-1 + yf G(xA, xA; x, y) G(a?i , A, 0; 0, x ) + ~ 4?1 " Ll±ij G{xX, xA; x, y) iJ(0; x ) 

+ 47f2( ~ 1 + y) G(xA, xA; x, y) H (1; x ) + 12 (-1 + yf G(xX, xA; x, y) H(0, 0, 0; y) 

- 8 (-1 + yf G(xA, xA; x, y) if (0, 1, 0; x ) - 8 (-1 + yf G(xA, xA; x, y) if (1, 0, 0; x ) 
+ 8 (-1 + yf G(xA, xA; x, y) if (1, 1, 0; x ) + 24 G(A, 0, A; x, y) if (0, 0; x ) 

- 4 G(A, 0, A; x, y) if (0, 0; y) + 8 G(A, 0, A; x, y) (1, 0; x ) 

+ (4 - 4 y) G(A, 0, xA; x, y) if (0, 0; y) + 8 (-1 + y) G(A, 0, xA; x, y) if (1, 0; y) 

+ 48 G(A, A, 0; x, y) H(0, 0; x ) - 20 G(A, A, 0; x, y) H(0, 0; y) 

+ 8 G(A, A, 0; x, y) if (1, 0; x ) + 16 G(A, x , A; x, y) if (0, 0; x ) 

+ 4 G(A, x , A; x, y) if (0, 0; y) + 16 G(A, x , A; x, y) if (1, 0; x ) 

+ 4 (-1 + y) G(A, x , xA; x, y) if (0, 0; y) + 8 (-1 + y) G(A, x , xA; x, y) if (1, 0; y) 

+ 16 G(A, xi, A; x, y) if (0, 0; x ) + 4 G(A, xi , A; x, y) if (0, 0; y) 

+ 16 G(A, xi, A; x, y) if (1, 0; x ) + 4 (-1 + y) G(A, x 1; xA; x, y) if (0, 0; y) 
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+ 8 (-1 + y) G(A, x 1: xX; x, y) if (1, 0; y) + (8 - 8 y) G(A, xX, 0; x, y) if (0, 0; y) 
+ 8 (-1 + y) G(A, xA, 0; x, y) if (1, 0; y) - 16 G(x , A, A; x, y) if (0, 0; x ) 

- 4 G(x , A, A; x, y) if (0, 0; y) - 16 G(x , A, A; x, y) if (1, 0; x ) 

+ (4 - Ay) G(x , A, xA; x, y) if (0, 0; y) + (8 - 8 y) G(x , A, xA; x, y) if (1, 0; y) 
+ (4 - 4y) G(x , xX, A; x, y) if (0, 0; y) + (8 - 8 y) G(x , xX, A; x, y) if (1, 0; y) 

- 4 (-1 + y) 2 G(x , xA, xA; x, y) ff(0, 0; y) - 16 G(xi, A, A; x, y) if (0, 0; x ) 

- 4 G(x x , A, A; x, y) if (0, 0; y) - 16 G(x 1; A, A; x, y) fi(l, 0; x ) 

+ (4 - 4 y) G(xi, A, xA; x, y) if (0, 0; y) + (8 - 8 y) G(x 1; A, xA; x, y) if (1, 0; y) 
+ (4 - 4 y) G(xi, xX, A; x, y) if (0, 0; y) + (8 - 8 y) G(xi, xA, A; x, y) if (1, 0; y) 

- 4 (-1 + yf G(xx , xA, xA; x, y) fi(0, 0; y) + 8 (-1 + y) G(xA, 0, A; x, y) ff (0, 0; y) 
+ 16 (-1 + y) G(xA, 0, A; x, y) if (1, 0; y) + 8 (-1 + yf G(xA, 0, xA; x, y) if (0, 0; y) 
+ (8 - 8 y) G(xA, A, 0; x, y) if (0, 0; y) + 8 (-1 + y) G(xA, A, 0; x, y) if (1, 0; y) 

+ 4 (-1 + y) G(xA, x , A; x, y) if (0, 0; y) + 8 (-1 + y) G(xA, x , A; x, y) if (1, 0; y) 
+ 4 (-1 + yf G(xA, x , xA; x, y) if (0, 0; y) + 4 (-1 + y) G(xA, xi, A; x, y) if (0, 0; y) 
+ 8 (-1 + y) G(xA, X! , A; x, y) if (1, 0; y) + 4 (-1 + yf G(xA, an , xA; x, y) if (0, 0; y) 
+ 4 (-1 + yf G(xA, xA, 0; x, y) if (0, 0; y) - 8 (-1 + yf G(xA, xA, 0; x, y) if (1, 0; x ) 
+ 16 G(A, 0, 0, A; x, y) if (0; x ) - 8 G(A, 0, 0, A; x, y) ff (0; y) 
+ (8 - 8 y) G(A, 0, 0, xA; x, y) if (0; y) + 24 G(A, 0, A, 0; x, y) if (0; x ) 

- 4 G(A, 0, A, 0; x, y) if (0; y) + 8 G(A, 0, x , A; x, y) if (0; x ) 

- 4 G(A, 0, x Q , A; x, y) if (0; y) + (4 - 4 y) G(A, 0, x , xA; x, y) if (0; y) 
+ 8 G(A, 0, x u A; x, y) if (0; x ) - 4 G(A, 0, an, A; x, y) if (0; y) 

+ (4 - 4 y) G(A, 0, x u xX; x, y) if (0; y) + (4 - 4 y) G(A, 0, xA, 0; x, y) if (0; y) 
+ 48 G(A, A, 0, 0; x, y) if (0; x ) - 8 G(A, A, 0, 0; x, y) if (0; y) 

- 72 G(A, A, A, A; x, y) if (0; x ) + 36 G(A, A, A, A; x, y) if (0; y) 

+ 36 (-1 + y) G(A, A, A, xA; x, y) if (0; y) + 8 G(A, x , 0, A; x, y) if (0; x ) 

- 4 G(A, x , 0, A; x, y) if (0; y) + (4 - 4 y) G(A, x , 0, xA; x, y) if (0; y) 
+ 16 G(A, x , A, 0; x, y) if (0; x ) - 8 G(A, x , A, 0; x, y) if (0; y) 

+ (4 - 4 y) G(A, x Q , xA, 0; x, y) if (0; y) + 8 G(A, x u 0, A; x, y) if (0; x ) 

- 4 G(A, xi, 0, A; x, y) if (0; y) + (4 - 4 y) G(A, x x , 0, xA; x, y) if (0; y) 
+ 16 G(A, x u A, 0; x, y) if (0; x ) - 8 G(A, x u A, 0; x, y) if (0; y) 

+ (4 - 4 y) G(A, xx , xA, 0; x, y) if (0; y) + 24 (-1 + y) 2 G(A, xA, xA, A; x, y) if (0; x ) 

- 12 (-1 + y) 2 G(A, xA, xA, A; x, y) if (0; y) - 12 (-1 + y) 3 G(A, xA, xA, xA; x, y) if (0; y) 

- 8 G(x , A, 0, A; x, y) if (0; x ) + 4 G(x , A, 0, A; x, y) if (0; y) 

+ 4 (-1 + y) G(x , A, 0, xA; x, y) if (0; y) - 16 G(x , A, A, 0; x, y) if (0; x ) 
+ 8 G(x , A, A, 0; x, y) ff (0; y) + 4 (-1 + y) G(x , A, xA, 0; x, y) ff (0; y) 
+ 4 (-1 + y) G(x ,xA, A,0;x,y)if(0;y) - 8 G{x u X, 0, A; x, y) ff (0; x ) 
+ 4 G(xi , A, 0, A; x, y) if (0; y) + 4 (-1 + y) G(xi , A, 0, xA; x, y) ff (0; y) 

- 16 G(xu A, A, 0; x, y) ff (0; x ) + 8 G(xi , A, A, 0; x, y) ff (0; y) 

+ 4 (-1 + y) G(xi, A, xA, 0; x, y) ff (0; y) + 4 (-1 + y) G(xi, xA, A, 0; x, y) ff (0; y) 
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+ (8 - 8 y) G(xX, 0, A, 0; x, y) ff (0; y) + 24 (-1 + yf G(x\, A, xA, A; x, y) ff(0; x ) 

- 12 (-1 + yf G(xX, A, xX, A; x, y) ff (0; y) - 12 (-1 + y) 3 G(xA, A, xX, xA; x, y) ff (0; y) 
+ (4 - Ay) G(xA, x , A, 0; x, y) ff(0; y) + (4 - Ay) G(xX, an, A, 0; x, y) if (0; y) 

+ 24 (-1 + yf G(xX, xX, A, A; x, y) ff (0; x ) - 12 (-1 + yf G(xX, xX, A, A; x, y) ff (0; y) 

- 12 (-1 + yf G(xX, xX, A, xA; x, y) ff (0; y) + 4 AG(0; 0, x ) G(A, 0; 0, x ) G(A, A; x, y) 

- 4 (-1 + yf AG(0; 0, x ) G(A, 0; 0, x ) G(xA, xA; x, y) + 4 AG(0; 0, x ) G(A, A; x, y) ff (0, 0; x ) 

- 4 AG(0; 0, x ) G(A, A; x, y) if (1, 0; x ) - 4 (-1 + yf AG(0; 0, x ) G(xA, xA; x, y) if (0, 0; x ) 
+ 4 (-1 + y) 2 AG(0; 0, x ) G(xA, xA; x, y) ff (1, 0; x ) + 4 AG(0; 0, x ) G(A, 0, A; x, y) ff (0; x ) 
+ 4 AG(0; 0, x ) G(A, A, 0; x, y) ff (0; x ) + 8 AG(0; 0, x ) G(A, x , A; x, y) if (0; x ) 

+ 8 AG(0; 0, x ) G(A, xi , A; x, y) ff (0; x ) - 8 AG(0; 0, x ) G(x , A, A; x, y) if (0; x ) 

- 8 AG(0; 0, x ) G(xi, A, A; x, y) ff (0; x ) - 4 (-1 + y f AG(0; 0, x ) G(xA, xA, 0; x, y) ff (0; x ) 
+ 8 AG(0; 0, y) G(A, A; x, y) if (0, 0; y) + 12 (-1 + y) AG(0; 0, y) G(A, xA; x, y) ff(0, 0; y) 

+ (8 - 8 y) AG(0; 0, y) G(A, xA; x, y) if (1, 0; y) + 12 (-1 + y) AG(0; 0, y) G(xA, A; x, y) if (0, 0; y) 
+ (8 - 8 y) AG(0; 0, y) G(xA, A; x, y) if (1, 0; y) + 4 (-1 + y) AG(0; 0, y) G(A, 0, xA; x, y) ff (0; y) 
+ 4 + AG(0;0,y)G(A,x ,xA;x,y)ff(0;y) + 4 (-1 + y) AG(0; 0, y) G(A, an, xA; x, y) ff(0; y) 
+ 4 (-1 + y) AG(0; 0, y) G(A, xA, 0; x, y) ff (0; y) + (4 - 4y) AG(0; 0, y) G(x , A, xA; x, y) ff (0; y) 
+ (4 - 4 y) AG(0; 0, y) G(x , xA, A; x, y) ff (0; y) + (4 - 4 y) AG(0; 0, y) G(an , A, xA; x, y) ff (0; y) 
+ (4 - 4 y) AG(0; 0, y) G(x u xX, A; x, y) ff (0; y) + 8 (-1 + y) AG(0; 0, y) G(xA, 0, A; x, y) ff (0; y) 
+ 4 (-1 + y) AG(0; 0, y) G(xA, A, 0; x, y) ff (0; y) + 4 (-1 + y) AG(0; 0, y) G(xA, x , A; x, y) ff (0; y) 
+ 4 (-1 + y) AG(0; 0, y) G(xA, an, A; x, y) ff (0; y) + 2 AG(0, 0; 0, x ) G(A, A; x, y) ff (0; x ) 

- 2 (-1 + yf AG(0, 0; 0, x ) G(xA, xA; x, y) ff (0; x ) + (4 - 4 y) AG(0, 0; 0, y) G(A, xA; x, y) ff (0; y) 
+ (4 - 4 y) AG(0, 0; 0, y) G(xA, A; x, y) ff (0; y) + 6 AG(0, A; 0, x ) G(A, A; x, y) ff (0; x ) 

- 12 AG(0, A; 0, x ) G(A, A; x, y) ff (0; y) - 6 (-1 + y f AG(0, A; 0, x ) G(xA, xA; x, y) ff (0; x ) 
+ 12 (-1 + yf AG(0, A; 0, x ) G(xA, xA; x, y) ff (0; y) + 8 AG(0, A; 0, y) G(A, A; x, y) ff (0; y) 

- 16 (-1 + y) AG(0, A; 0, y) G(A, xA; x, y) ff (0; x ) + 8 (-1 + y) AG(0, A; 0, y) G(A, xA; x, y) ff (0; y) 

- 16 (-1 + y) AG(0, A; 0, y) G(xA, A; x, y) ff (0; x ) + 8 (-1 + y) AG(0, A; 0, y) G(xA, A; x, y) ff (0; y) 
+ 6 (-2 Vy + y) AG(0, xA; 0, x ) G(A, A; x, y) ff (0; x ) 

- 6 (-2 + Vy) (-1 + y) 2 VV A G(0, xA; 0, x ) G(xA, xA; x, y) ff (0; x ) 

+ 4 (-1 + yf AG(0, xA; 0, y) G(A, xA; x, y) ff (0; y) + 4 (-1 + yf AG(0, xA; 0, y) G(xA, A; x, y) ff(0; y) 
+ 4 G(A; 0, x ) G(A, 0; 0, x ) G(A, A; x, y) - 4 (-1 + y) 2 G(A; 0, x ) G(A, 0; 0, x ) G(xA, xA; x, y) 
+ 4 G(A; 0, x ) G(A, A; x, y) ff (0, 0; x ) - 16 G(A; 0, x ) G(A, A; x, y) ff (0, 0; y) 

- 4 G(A; 0, x ) G(A, A; x, y) ff (1, 0; x ) - 4 (-1 + y) 2 G(A; 0, x ) G(xA, xA; x, y) ff (0, 0; x ) 

+ 16 (-1 + yf G(A; 0, x ) G(xA, xA; x, y) ff (0, 0; y) + 4 (-1 + y) 2 G(A; 0, x ) G(xA, xA; x, y) ff (1, 0; x ) 
+ 4 G(A; 0, x ) G(A, 0, A; x, y) ff (0; x ) - 8 G(A; 0, x ) G(A, 0, A; x, y) ff (0; y) 
+ 4 G(A; 0, x ) G(A, A, 0; x, y) ff (0; x ) - 8 G(A; 0, x ) G(A, A, 0; x, y) ff (0; y) 
+ 8 G(A; 0, x ) G(A, x , A; x, y) ff (0; x ) - 16 G(A; 0, x ) G(A, x , A; x, y) ff (0; y) 
+ 8 G(A; 0, x Q ) G(A, an, A; x, y) ff (0; x ) - 16 G(A; 0, x ) G(A, an, A; x, y) ff (0; y) 

- 8 G(A; 0, x ) G(x , A, A; x, y) ff (0; x ) + 16 G(A; 0, x ) G(x , A, A; x, y) ff (0; y) 
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- 8 G(A; 0, x ) G{x x , A, A; x, y) ff (0; x ) + 16 G(A; 0, x ) G{ Xl , A, A; x, y) ff (0; y) 

- 4 (-1 + yf G(A; 0, x ) G(x\, xX, 0; x, y) ff (0; x Q ) + 8 (-1 + yf G(A; 0, x ) G(xA, xX, 0; x, y) ff (0; y) 
+ 8 G(A; 0, y) G(A, A; z, y) ff (0, 0; y) + 16 G(A; 0, y) G(A, A; x, y) ff (1, 0; y) 

- 32 (-1 + y) G(A; 0, y) G(A, xA; x, y) ff (0, 0; x ) + 12 (-1 + y) G(A; 0, y) G(A, xA; x, y) ff (0, 0; y) 

- 32 (-1 + y) G(A; 0, y) G(xA, A; x, y) ff (0, 0; x ) + 12 (-1 + y) G(A; 0, y) G(xA, A; x, y) ff (0, 0; y) 
+ (8 - 8 y) G(A; 0, y) G(A, 0, xA; x, y) if (0; x ) + 4 (-1 + y) G(A; 0, y) G(A, 0, xA; x, y) if (0; y) 

+ (8 - 8 y) G(A; 0, y) G(A, x , xA; x, y) if (0; x ) + 4 (-1 + y) G(A; 0, y) G(A, x , xA; x, y) if (0; y) 
+ (8 - 8 y) G(A; 0, y) G(A, x 1; xA; x, y) ff(0; x ) + 4 (-1 + y) G(A; 0, y) G(A, x 1; xA; x, y) ff(0; y) 
+ (8 - 8 y) G(A; 0, y) G(A, xA, 0; x, y) if (0; x ) + 4 (-1 + y) G(A; 0, y) G(A, xA, 0; x, y) ff (0; y) 
+ 8 (-1 + y) G(A; 0, y) G(x , A, xA; x, y) ff (0; x ) + (4 - 4 y) G(A; 0, y) G(x , A, xA; x, y) ff (0; y) 
+ 8 (-1 + y) G(A; 0, y) G(x , xA, A; x, y) ff(0; x ) + (4 - 4 y) G(A; 0, y) G(x , xA, A; x, y) ff(0; y) 
+ 8 (-1 + y) G(A; 0, y) G(xi, A, xA; x, y) ff (0; x ) + (4 - 4y) G(A; 0, y) G(xi, A, xA; x, y) ff (0; y) 
+ 8 (-1 + y) G(A; 0, y) G(x u xX, A; x, y) ff(0; x ) + (4 - 4y) G(A; 0, y) G(x u xX, A; x, y) ff(0; y) 

- 16 (-1 + y) G(A; 0, y) G(xA, 0, A; x, y) ff (0; x ) + 8 (-1 + y) G(A; 0, y) G(xA, 0, A; x, y) ff (0; y) 
+ (8 - 8 y) G(A; 0, y) G(xA, A, 0; x, y) ff (0; x ) + 4 (-1 + y) G(A; 0, y) G(xA, A, 0; x, y) ff (0; y) 

+ (8 - 8 y) G(A; 0, y) G(xA, x , A; x, y) ff (0; x ) + 4 (-1 + y) G(A; 0, y) G(xA, x , A; x, y) ff (0; y) 
+ (8 - 8y) G(A; 0, y) G(xA, xi, A; x, y) ff(0; x ) + 4 (-1 + y) G(A; 0, y) G(xA, x u A; x, y) ff(0; y) 
+ 2 G(A, 0; 0, x ) G(A, A; x, y) if (0; x ) - 16 G(A, 0; 0, x ) G(A, A; x, y) if (0; y) 

- 2 (-1 + y) 2 G(A, 0; 0, x ) G(xA, xA; x, y) ff (0; x ) + 16 (-1 + yf G(A, 0; 0, x ) G(xA, xA; x, y) ff (0; y) 

- 24 (-1 + y) G(A, 0; 0, y) G(A, xA; x, y) ff (0; x ) + (4 - 4 y) G(A, 0; 0, y) G(A, xA; x, y) ff (0; y) 

- 24 (-1 + y) G(A, 0; 0, y) G(xA, A; x, y) ff (0; x ) + (4 - 4 y) G(A, 0; 0, y) G(xA, A; x, y) H(0; y) 
+ 6 G(A, A; 0, x ) G(A, A; x, y) if (0; x ) - 12 G(A, A; 0, x ) G(A, A; x, y) ii(0; y) 

- 6 (-1 + y) 2 G(A, A; 0, x ) G(xA, xA; x, y) if (0; x ) + 12 (-1 + y) 2 G(A, A; 0, x ) G(xA, xA; x, y) if (0; y) 
+ 16 G(A, A; 0, y) G(A, A; x, y) if (0; x ) - 8 G(A, A; 0, y) G(A, A; x, y) if (0; y) 

+ 6 (-2Vy + ?/) G(A,A;x,y)G(A,xA;0,x o )if(0;x o ) 

- 16 (-1 + y) G(A, A; x, y) G(A, xA; 0, y) if (0; y) 

+ 4 G(A, A; x, y) G(x , A; 0, x ) if (0; x ) - 8 G(A, A; x, y) G(x , A; 0, x ) if (0; y) 
+ 4 (-2Vy + y) G(A,A;x,y)G(x o ,xA;0,x o )if(0;xo) 

+ 4 G(A, A; x, y) G(xi , A; 0, x ) if (0; x ) - 8 G(A, A; x, y) G(xi , A; 0, x ) if (0; y) 
+ 4 (-2 + y) G(A, A; x, y) G(xi , xA; 0, x ) if (0; x ) 

- 4 G(A, A; x, y) if (0; x ) if (0, 0; y) - 4 G(A, A; x, y) if (0; y) if (0, 0; x ) 

- 6 (-2 + y/y) (-1 + y) 2 x /yG(A, xA; 0, x ) G(xA, xA; x, y) if (0; x ) 

+ 8 (-1 + y) 2 G(A, xA; 0, y) G(A, xA; x, y) if (0; x ) - 4 (-1 + y) 2 G(A, xA; 0, y) G(A, xA; x, y) if (0; y) 
+ 8 (-1 + y) 2 G(A, xA; 0, y) G(xA, A; x, y) if (0; x ) - 4 (-1 + y) 2 G(A, xA; 0, y) G(xA, A; x, y) if (0; y) 

- 16 (-1 + y) G(A,xA;x,y)G(x o ,A;0,y)if(0;x o ) +8 (-1 + y) G(A, xA; x, y) G(x , A; 0, y) ff (0; y) 
+ 8 (-1 + y) 2 G(A, xA; x, y) G(x , xA; 0, y) if (0; y) - 16 (-1 + y) G(A, xA; x, y) G(xx , A; 0, y) if (0; x ) 
+ 8 (-1 + y) G(A, xA; x, y) G(xi , A; 0, y) ff (0; y) + 8 (-1 + y) 2 G(A, xA; x, y) G(xi , xA; 0, y) ff (0; y) 
+ 4 (-1 + y) G(A, xA; x, y) ff (0; x ) ff (0, 0; y) + 4 (-1 + y) G(A, xA; x, y) ff (0; y) ff (0, 0; x ) 

- 4 (-1 + y) 2 G(x , A; 0, x ) G(xA, xA; x, y) ff (0; x ) + 8 (-1 + yf G(x , A; 0, x ) G(xA, xA; x, y) ff (0; y) 
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-16 (-1 + y) G(x ,A;0,y)G(xA,A;x,y)ff(0;x )+8 (-1 + y) G(x , A; 0, y) G(x\, A; x, y) JJ(0; y) 

- 4 (-2 + Vy) (-1 + y) 2 G(x , xX; 0, ar ) G(xA, xA; x, y) Jf (0; x ) 

+ 8 (-1 + yf G(x ,xX; 0, y) G(xA, A; x, y) H(0; y) - 4 (-1 + y) 2 G(x : , A; 0, x ) G(xX, xA; x, y) H(0; x ) 
+ 8 (-1 + yf G{ Xl , A; 0, x ) G(xX, xA; x, y) H(0; y) - 16 (-1 + y) G{x u A; 0, y) G(xA, A; x, y) (0; x ) 
+ 8 (-1 + y) G( Xl ,X;0,y)G(xX,X;x,y) H(0;y) 

- 4 (-2 + ^y) (-1 + vf y/yG(xi,x\\ 0, ar ) G(xX, xA; x, y) H(0; x ) 

+ 8 (-1 + yf G( Xl , xA; 0, y) G(xX, A; x, y) H(0; y) + 4 (-1 + y) G(xX, A; x, y) ff(0; x ) H(0, 0; y) 
+ 4 (-1 + y) G(xA, A; x, y) H(0; y) H(0, 0; x ) - 4 (-1 + yf G(xX, xA; x, y) (0; x ) H (0, 0; y) 

- 4 (-1 + yf G(xX, xA; x, y) if (0; y) if (0, 0; x ) - 4 G(A, A, 0; x, y) H(0; x ) if (0; y) 

+ 4 (-1 + y) G(A, xA, 0; x, y) H(0; x ) (0; y) + 4 (-1 + y) G(xA, A, 0; x, y) if (0; x ) H(0; y) 

- 4 (-1 + y) 2 G(xA, xA, 0; x, y) H(0; x ) H (0; y). 

As with the other MI, we have checked that the leading (finite) contribution agrees with the results 
of Rcf. [54]. 

7. Summary 

In this paper, we have provided series expansions in the dimensional regularisation parameter e for 
all two- loop Master Integrals with three external off-shell legs and all internal lines being masslcss. 
The results are presented in terms of an extended basis of 2-dimensional harmonic polylogarithms. 
The novel feature is that this basis includes quadratic forms - that matches on to the allowed phase 
space boundary for the 1 — > 2 decay. For each Master Integral, we have given sufficient terms in 
the e-expansion to describe two-loop vertex corrections for physical processes. 

The MI presented here are ingredients for a variety of interesting two-loop processes such as 
the QCD corrections to H — ► V*V* decay in the heavy top quark limit and the QCD corrections 
to the fully off-shell triple gluon (and quark-gluon) vertices. 

The MI also form a staging post for the study of massless two-loop 2 — > 2 scattering amplitudes 
with two off-shell legs. These processes include the NNLO QCD corrections to qq — ► V*V* (where 
V = W,Z) and the NLO corrections to gg — > V*V*. Altogether there are 11 planar box and 3 
non-planar box master topologies which remain to be studied. 
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